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1. Abstract 

Thousands of extrasolar planets have been discovered, and it is clear that the galactic 
planetary census draws on a diversity greatly exceeding that exhibited by the solar system’s 
planets. We review signihcant landmarks in the chronology of extrasolar planet detection, 
and we give an overview of the varied observational techniques that are brought to bear. We 
then discuss the properties of the planetary distribution that is currently known, using the 
mass-period diagram as a guide to delineating hot Jupiters, eccentric giant planets, and a 
third, highly populous, category that we term “ungiants”, planets having masses M < 30 
and orbital periods P < 100 d. We then move to a discussion of the bulk compositions 
of the extrasolar planets, with particular attention given to the distribution of planetary 
densities. We discuss the long-standing problem of radius anomalies among giant planets, 
as well as issues posed by the unexpectedly large range in sizes observed for planets with 
mass somewhat greater than Earth’s. We discuss the use of transit observations to probe 
the atmospheres of extrasolar planets; various measurements taken during primary transit, 
secondary eclipse, and through the full orbital period, can give clues to the atmospheric 
compositions, structures and meteorologies. The extrasolar planet catalog, along with the 
details of our solar system and observations of star-forming regions and protoplanetary disks, 
provide a backdrop for a discussion of planet formation in which we review the elements of 
the favored pictures for how the terrestrial and giant planets were assembled. We conclude 
by listing several research questions that are relevant to the next ten years and beyond. 
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Introduction 


Until nearly the close of the twentieth centnry, the search for planets beyond the Solar 
System was a franght and generally less-than-fnlly respectable endeavor. For generations, 
the scientihc commnnity was snbjected to spnrious annonncements of discovery of planets 
orbiting nearby stars. Famous examples include erroneous claims, starting in 1855 of a 


planetary mass object in the 70 Ophiuchi system (e.g. 

Jacob 

i.s55 

See 

1896 

), and that one. 

or possibly two, Jupiter-like planets were harbored by Barnard’s 

Star, the red dwarf that 


is our Sun’s nearest stellar neighbor aside from the a Centauri triple-star system (van de 


Kamp 1963, 1969). 


If we adopt the frequently-used definition that requires the planetary mass, Mp, to lie 
below the 13 Mj deuterium-burning limit, where Mj denotes the mass of Jupiter, then it 


might be appropriate to credit the first bona-hde detection of an extrasolar planet to Latham 


et ah (1989), who used Doppler radial velocity measurements to sense an object with a mini¬ 


mum mass of 11 Jupiter masses on an eccentric 83.9-day orbit around HD 114762, a relatively 
nearby, solar-type star. With hindsight, this object belongs to the class of giant planets that 
tend to be more massive than Jupiter, which often have substantial orbital eccentricities, and 
which are found in orbit around approximately 5% of local solar-type dwarfs (Butler et ah 


2006). At the time of HD 114762 b’s discovery, however, as a consequence of the abysmal 


track record of extrasolar planet claims, Latham et ah’s announcement did not generate 
signihcant excitement within the wider community. Indeed, the profound difference between 
HD 114762 b’s mass and orbit from those of our Solar System’s planets led Latham et ah to 
state “This leads to the suggestion that the companion is probably a brown dwarf, and 
may even be a giant planet.” Given that HD 114762 b’s product of its mass and sine of the 
inclination angle between its orbital plane, Mpsini, is so near this physical mass boundary 
and sini < 1, Latham et ah’s statement can still be considered as a proper assessment of 
our current knowledge of this object a quarter century later! 

Soon thereafter, astronomers developed a sudden and abiding interest in extrasolar 


planets. The catalyst was the pulse-timing detection, by Wolszczan & Frail (1992), of two 


terrestrial-mass planets, with Mpsini = 3.4 M 0 and 2.8 M 0 , and with Mercury-like periods 
of 66.5 and 98.2 days, orbiting the millisecond pulsar PSR B1257-I-12. Pulsars, which are 
rapidly rotating, strongly magnetized neutron stars, emit radio waves that appear as periodic 
pulses to an observer on Earth. Although a pulsar emits an exquisitely periodic signal, the 
times at which the pulses reach the receiver are not equally spaced if the distance between the 
pulsar and the telescope varies in a nonlinear fashion, and if periodic variations are present in 
data that has been properly reduced, they may indicate the presence of companions orbiting 
the pulsar. In the case of PSR B1257-I-12, measurements are so accurate that deviations 
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from keplerian orbits of the two planets mentioned above cansed by their mutual interactions 
have been measured and used to calculate their actual masses, Mp = 4.3 ± 0.2 M 0 and 3.9 
± 0.2 M®, implying that the system is inclined by ~ 50° from the line of sight (Konacki & 


Wolszcz^|2003 ). A handful of other pulsar planets have been found, including a lunar-mass 
third planet orbiting PSR B1257-I-12, and a long-period ~ 2.5Mj planet in a circumbinary 


orbit about a pulsar-white dwarf binary system in the M4 globular cluster (Sigurdsson et ah 


2003). 


The pulsar planets were followed in 1995 by Michel Mayor and Didier Queloz’s discovery 
of 51 Pegasi bQ a gas giant planet with a minimum mass roughly half that of Jupiter, and a 
startlingly short 4.2308-day orbital period about a metal-rich, but otherwise ordinary G-type 
main sequence star. Thereafter, the pace of discovery ballooned, along with our understand¬ 
ing of both the planetary census and of the physical characteristics of alien worlds. The 
study of exoplanet physical characteristics has largely been made possible by the detection 
of transiting planets, which partially eclipse the light of their star once per orbit, the first of 


which was discovered in 1999 in orbit around HD 209458 (Henry et al. 2000 Charbonneau 


et al. 2000). 


It is no exaggeration to say that a fully independent branch of astronomy has been 
created from scratch. The rapid pace of progress is reflected in Figure 1, which charts 
log^Q Mp sin z as a function of year of discovery. The lower envelope of the discovery masses 
in this semi-log plot is well described by a constant slope, indicative of a Moore’s-law behavior 
with a r ~ 2 year halving time of the lowest-known extrasolar planetary mass. 


^ “Official” exoplanetary names are a compound of the name of the parent star and a lower-case letter. 
The first planet detected in orbit around a given star is denoted 6, the second, c, etc. 
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Fig. 1.— Discovery date vs. estimated masses for extrasolar planets orbiting normal stars 
(not pnlsars) over time, with a logarithmic ?/-axis showing estimated masses (or for the cases 
of Doppler detections, Msini’s), and year of discovery rnnning along the x-axis. The lower 
envelope of this semi-log plot shows a clear Moore’s Law-like progression, and demonstrates 
that the held has collectively arrived (in the past 1-2 years) at the discovery of planets 
with ~ 1 M 0 . Transiting planets are shown in red, whereas planets detected only via 
Doppler radial velocity are shown in gray. Data are taken from www.exoplanets.org/csv- 
hles/exoplanets.csv (Han et ah 2014). 






































3. Planetary Detection Methods 


Given the inherent difficulties associated with the detection of extrasolar planets, it 
is remarkable that a very broad range of observational methods for characterization and 
detection have been either successfully employed, or are being seriously prepared for future 
use. As a corollary to the vast interstellar distances, all extrasolar planets are extremely faint; 
to a distant optical observer, Jupiter’s brightness is Lj ~ 1 x 10“® the solar value in the V 
band (although it rises to a contrast of Aj/L© ~ 2 x 10“"^ redward of 30 /rm). Most detection 
methods, therefore, are indirect, in that the planets are sensed through their influence on the 
host star or on the path of light emitted by a distant star. The range of distinct observation 
methods are sensitive to different classes of planets, and importantly, they provide us with 
complementary information about the planets they do hnd. Without question, a full quiver 
of techniques is required in order to provide an optimal characterization of the diversity of 
planets and of planetary systems. 

To date, two methods - photometric transit detection and Doppler radial velocity mea¬ 
surement - have generated most of the successful detections. The Doppler radial velocity 
(RV) technique hinges on the ability to monitor tiny time-dependent wavelength displace¬ 
ments among the profusion of atomic lines in a star’s spectrum. By measuring the net overall 
shift, X = SX/X, using a large number of individual spectral lines, the velocity at which the 
star is moving towards or away from an observer can be precisely determined. Astronomers 
then subtract the motion of the telescope relative to the barycenter of the Solar System and 
other known motions, revealing the radial motion of the target star that results from the tug 
of its own planets. Transit photometry identihes planet candidates by the periodic drops in 
observed stellar brightness that are caused by the planet obscuring a portion of the stellar 
disk once per orbit. Note that transit photometry can only detect planets whose orbits are 
viewed nearly edge-on. 

Exoplanet detection using either the radial velocity method or the transit photome¬ 
try method requires very precise measurements of light from the planet-hosting star. The 
techniques used to achieve such high precision are based on differential measurements of ob¬ 
servations taken at different times. Absolute accuracies, therefore, need not be particularly 
high, so the mean velocity (relative to the Solar System) of a star hosting RV-detected plan¬ 
ets is generally not known to an accuracy approaching 1 ms“^, nor is the spectrophotometric 
brightness of a star hosting a planet observed via transit photometry known to the, e.g., part 
per thousand level needed in relative brightness to detect Neptune-size planets. 

As with pulsar timing, radial velocity measurements yield the product of the planet’s 
mass (divided by that of the star, whose mass can usually be estimated fairly accurately 
from its spectral characteristics) and the sine of the inclination angle between the orbital 
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plane and the plane of the sky, Mpsinh The total half-amplitnde, K, of the stellar reflex 
velocity indnced by an orbiting planet is the measnred qnantity, and is related to Mp sin i 
throngh 

Mpsini 1 

P J (M* + MpfP Vl-eP 


K = 


( 1 ) 


indicating that the radial velocity techniqne is most sensitive to massive planets and to 
planets in short-period orbits. The RV techniqne allows estimation of the period, P, orbital 
eccentricity, e, the mean anomaly, M = 27r(T — Tperiastron)/.P, and orientation - generally 
expressed as the longitnde of periapse, oj, of the stellar orbit. Mnlti-parameter estimation 
techniqnes have reached a high degree of rehnement for the problem of determining best-£t 
orbital elements and nncertainty estimates. For reviews that cover the process of radial 
velocity modeling, see Ford (2006); Meschiari et ah ( 2009| ); Dawson & Fabrycky (2010); 
Dnmnsqne et ah (2011)]^ 


Leading research gronps are now rontinely achieving instrnmentally-based velocity er¬ 
rors of (j ~ 50 cm s“^ (representing a Doppler shift of two parts in 10®) on stable, bright, 
chromospherically inactive stars. For context, Jnpiter canses the Snn’s velocity to vary with 
a half-amplitnde Kj = 12.5 ms“^ and a period Pj = 11.86 years; Satnrn’s effect is the next 
largest, with an amplitnde of 2.8 ms“^ and a period, Ps, of almost 30 years, bnt Earth in- 
dnces a mere = 9 cms“^ flnctnation. To date, the planet with the lowest radial velocity 
amplitnde signal yet annonnced is Alpha Centanri B b, with K = 51 cms“\ P = 3.2357 d. 


and Msini = 1.1 M 0 (Dnmnsqne et al. 2012). Earth-like planets orbiting snn-like stars at 


astrobiologically interesting distances (where liqnid water conld be stable on their snrfaces) 
remain somewhat beyond the capabilities cnrrently achieved. 

Precise radial velocity measnrements reqnire the good statistics that come from aggre¬ 
gating individnal measnrements of a large nnmber of spectral lines. High precision measnre¬ 
ments thns cannot be achieved for hotter stars with temperatnres Teg > 6000 K (correspond¬ 
ing to spectral types O, B, A, and early F), as these stars have far fewer spectral featnres 
than do stars of solar type or cooler. Stellar rotation and intrinsic variability (inclnding star 
spots) also represent major sonrces of noise that confonnd radial velocity measnrements. 
The radial velocity snrveys, snch as the Geneva Planet Search (Mayor et al.||2009), the Cali¬ 


fornia Planet Search (Howard et al. 2011), and the Lick-Carnegie Planet Search (Vogt et al. 


2010), are therefore increasingly focnsing on very nearby, old, photometrically qniet late G, 


K, and M-type main seqnence dwarfs. The 451 stars that constitnte the HARPS Gnaranteed 


^We recommend the Systemic Console software package, available at http;//www.stefanom.org/ 
systemic/, as an intuitive framework for exploring radial velocity and transit timing data. 
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Time Observations (GTO) high precision sample ( Sousa et ah||2008 ) encompass many of the 
Sun-like stars at declination, 6 < 20° N that are currently receiving intense surveillance. 

Planets can also be detected via the wobble that they induce in their parent star’s 
motion projected onto the plane of the sky. This so-called astrometric technique is most 
sensitive to massive planets orbiting stars that are relatively close to Earth. Because the 
star’s motion is detectable in two dimensions, the planet’s actual mass, rather than just 
the combination Msinz, can potentially be measured. Planets on more distant orbits are 
ultimately easier to detect using astrometry because the amplitude of the star’s motion 
is larger, but finding these planets requires a longer timeline of observations due to their 


greater orbital periods. To date, despite decades of optimistic expectations, see, e.g.. Black 


& Scargle (1982), the astrometric method has yet to unambiguously detect an exoplanet. 


and the most relevant discoveries are of objects such as the 28Mj companion to nearby LI.5 
dwarf DENIS-P J082303.1-491201 ( Sahlmann et al.|2013 ) that lies just above the upper mass 
limit for planets. The drought is likely to end soon, however. Simulations by|Casertano et 


(2008), indicated that ESA’s Gaia Mission, launched on December 19th, 2013, and designed 


to achieve 7 /ras precision on E = 10 stars, should be capable of detecting several thousand 
giant planets out to 3-4 AU from stars lying within 200 pc of the Sun. A more recent analysis 


by Perryman et al. (2014), which uses an updated galactic model and planet occurrence rates. 


and which folds in Gaia’s actual on-sky astrometric performance, suggests that the mission 
should hnd 21, 000 ± 6, 000 long-period planets with IMj < Mp < 15Mj during its nominal 
5-year duration, with several times more planets expected to be observed in the event that 
the mission is extended. 


The Doppler and astrometric techniques measure a planet’s gravitational pull on its 
star, and are thus sensitive to the planet’s mass. By contrast, transit photometry detects 
the amount of starlight that a planet can obscure, and thus yields an estimate of the planet’s 
size. If the line of sight from the parent star to the Earth lies in or near the orbital plane 
of an extrasolar planet, that planet passes in front of the disk of its star once each orbit as 
viewed from Earth, and periodically blocks a small fraction of the star’s light. Sufficiently 
precise photometric time-series measurements of the star’s brightness are used to reveal tran¬ 
sits, which can be distinguished from rotationally-modulated star spots and intrinsic stellar 
variability by their periodicity, their effectively square-well shapes, and their relative spec¬ 
tral neutrality. Although geometrical considerations limit the fraction of planets detectable 
by this technique (planetary orbital angular momentum vectors are aligned apparently ran¬ 
domly with respect to the galactic plane) thousands of stars can be surveyed within the 


^The angular momentum of the Solar System ultimately derives from the turbulent motions that existed 
within the Sun’s parent giant molecular cloud, and is not a direct consequence of the differential rotation of 
















field of view of one telescope, so transit photometry is potentially qnite efficient. Indeed, 
transit photometry now generates a clear majority of new planetary candidates, and gronnd- 
based snrveys snch as the HAT Net ( Bakos et al.||2012 ) and SnperWASP ( Smith et al.||2013 ), 
are generating dozens of conhrmed planets per year. 

The preliminary version of the Kepler Q12 planet candidate catalog, consisting of sig¬ 
nals identihed nsing the hrst 3 years of Kepler data, contains 3,538 candidate planets, at 
least ~ 90% of which are expected to represent true exoplanets (Morton & Johnson 2011 


Casertano et al. 

2008 

Santerne et al. 

2013 

Lissauer et al. 

2014b 


Almost half of these 

objects are associated with target stars having two or more planet candidates, implying 


that conhgurations of tightly spaced, highly co-planar planets are ubiquitous (Lissauer et ah 


2011b). More than 90% of these candidates have estimated radii Rp < Q R 0 and almost 
half have radii Rp < 2 R®, despite observational biases against detecting smaller planets]^ 
In many of Kepler’s multiple-transiting planet systems, the transit events exhibit deviations 
from strict periodicity. These so-called transit timing variations (as well as transit duration 
variations that are, much more rarely, also observed) can be attributed to planet-planet in¬ 
teractions within the systems, and have, in many cases, established to high conhdence that 
a particular candidate system is the genuine item rather than a false positive. Scores of 
Kepler candidates, furthermore, have been conhrmed with Doppler follow-up measurements, 
e.g., Cochran et al. ( |2011 ), Marcy et ah (2014), and hundreds of candidates have been vali¬ 
dated to be planets at very high conhdence through statistical arguments (e.g. Torres et al. 


2011 Rowe et al. 2014). Continued follow-up observations of Kepler Objects of Interest, as 


well as new photometric data from the K2 Mission (Howell et al. 2014) will keep the Kepler 
project productive for years to come. 

Reaching further into the box of astrophysical tricks, one hnds the microlensing method, 
which has been successfully employed to investigate the distribution of faint stellar and sub- 
stellar mass bodies within our galaxy ( Alcock et al.||1997| ), and which has now produced pub¬ 
lished detections of 18 extrasolar planets ( Kains et al.||2013 ). The microlensing phenomenon 
stems from the general relativistic bending of the light from a distant star (the source) by 
a massive object (the lens) passing between the source and the observer. Such lensing can 
cause the source to appear to gradually brighten by a factor of a few, or for special conhgu¬ 
rations leading to high magnihcation events factors exceeding 100, over a time scale of weeks 
to months. If the lensing star has planetary companions, then these far less massive bodies 


the galactic disk. 

^For an up-to-date list of Kepler candidates, see http://exoplcinetarchive.ipac.caltech.edu/ 
cgi-bin/ExoTables . 
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can produce brief enhancements in the brightness - generally of a few hours in duration - 
provided the line of sight passes close to the planet. Under favorable circumstances, planets 
as small as Earth can be detected. The properties of individual planets, however, and usu¬ 
ally even of the stars that they orbit, can only be estimated in a statistical sense because of 


the many parameters that influence a microlensing light curve see, e.g., Peale (1997); Gaudi 


(2012). The microlensing technique is comparatively well-suited to the detection of planets 
with long periods, and can readily probe Jovian, or even super-Earth mass planets in orbits 
similar to that of Jupiter itself. Due to the prevalence of M-type main-sequence stars in the 
stellar mass function, the microlensing technique has found a number of examples of giant 
planets and Neptune-mass planets beyond the locations of the so-called nebular snow lines 


for these stars (Bennett et ah 2010). Sumi et ah (2011) have extensively analyzed several 


years of microlensing survey data from the Galactic Bulge region, and report the discovery 
of a population of “isolated” planetary mass objects with M ~ Mj, which appear to occur 
with somewhat higher frequency than main-sequence stars. These planetary mass objects 
he either in very distant orbits {d > 100 AU) from undetected parent stars, or they are 
gravitationally unbound and thus floating freely through the galaxy. The recent discovery 
of WISE J085510.83071442.5 (Luhman 2014), a free-floating 3-10 Mj object lying a mere 2 
parsecs from the Sun, lends credence to the idea that free-floating planets could be common. 
Gertainly, the potential presence of hundreds of billions of such drifters, should these worlds 
withstand statistical scrutiny, would have important consequences for the theory of planet 
formation, and would serve as a dehnitive endorsement of the microlensing technique. 


All of the foregoing are indirect methods to locate extrasolar planets, but what about 
imaging planets directly? While observational efforts to directly detect extrasolar planets 
can be traced back to Ghristiaan Huygens in the 1600’^ only very recently have the ex¬ 
treme planet-star brightness contrasts and angular separations been successfully challenged. 
The reflected starlight from planets with orbits and sizes like those in our Solar System 
is well over twenty visual magnitudes (a factor > 10 ®) fainter than the stellar brightness. 
Diffraction of light by telescope optics and atmospheric variability compound the difficulty 
of direct detection of extrasolar planets. Technological advances, most critically in the area 
of adaptive optics on large ground-based telescopes and coronographic imaging from space, 
are enabling direct detection of relatively massive, young, and in most cases self-luminous, 
planets orbiting at substantial distances a > 10 AU from young, nearby stars with masses 


M* > Mq. These detections leverage planet-star brightness contrasts of order Lp/L* < 10^ 
and angular separations of order a ~ 1 arc sec. A high prohle example is illustrated by Beta 
Pictoris b (Lagrange et ah 2010), with mass Mp = 8^2 Mj at separation a = 8 IQ 4 AU from 


^Huygens, C. 1698 Cosmotheoros, Adriaan Moetjens: The Hague, Netherlands. 





















10 


Beta Pictoris, a young A-type star that has long been known to harbor an infrared excess 
and an extensive debris disk (Smith & Terrile 1984). Another interesting and somewhat 


puzzling instance is provided by Fomalhaut b ( Kalas et al.||2008 ). This curious object, which 
has been imaged with HST, lies just beyond the inner radius of the dust ring surrounding 
the young, nearby A-type star Fomalhaut. The imaged light associated with the planetary 
candidate’s position is dust-scattered starlight, associated, perhaps, with a circumplanetary 
disk (Galicher et ah 2013; Kenyon et ah 2014). 


The most startling and remarkable direct-imaging detection is that of the system of four 
companions to the young (~ 30Myr) A-type star HR 8799 ( Marois et al.||2008 2010). These 
planets all have masses in the vicinity Mp ~ 5 —7 Mj, with semi-major axes lying at 14.5, 24, 
38, and 68 AU, and call to mind a version of our own outer Solar System that has been scaled 
up by a factor of several in both mass and astrocentric distance. The presence of HR 8799’s 
planets at their observed distances presents interesting challenges to both the core-accretion 


(Lissauer et al. 2009) and the gravitational instability (Durisen et ah 2007) models of giant 


planet formation. The proximity of the system configuration to the dynamical stability 
boundary is also intriguing, and implies either that the system will suffer catastrophe within 
a time that is much less than the star’s main sequence life, or that the planets are dynamically 


protected by a resonant mechanism (Fabrycky & Murray-Clay 2010). The pace of detection 


by direct imaging is proceeding at a rapid clip, with additional examples being those of HD 


95086b (Rameau et al. 2013), 2MASS J01033563-5515561(AB)b (Rameau et al. 2013), Ksi 


And b (Carson et al. 2013), and GJ504b (Kuzuhara et al. 2013), which tend to be massive 


planets M > 4Mjup orbiting at i? > 40 AU from young (or relatively young) stars. 

Direct imaging is poised for major advances with the commissioning of several new 
projects that utilize 8 to 10-meter class telescopes. Among these are the Gemini Planet 
Imager (GPI), which experienced first light at the Gemini South Telescope in late 2013 
(Macintosh et al. 2014). The GPI adaptive optics system is designed to detect planetary 
companions at separations of 6=0.2 -1.0 arc seconds, and brightness contrasts of ~ 10“^ in 
1-2 hour exposures. This observational strategy is sensitive to Jupiter-mass planets within 
the hrst several million years of their formation, and will enable direct imaging of a variety of 
planets with ages under a billion years. The GPI system, furthermore, uses an integral field 
spectrometer to obtain i? ~ 40 spectra that will give important clues to overall atmospheric 
compositions and states. Already, spectra of directly imaged exoplanets are giving physical 


information, see, e.g. Konopacky et al. (2013), who report on the use of the Keck Telescope’s 


OSIRIS integral field spectrograph to detect GO and H2O absorption in HR 8799c. The 


SPHERE adaptive optics instrument on the VLT (Beuzit et al. 2008) draws on optical and 


near-IR imaging, low-resolution spectroscopy, and polarimetry in the service of finding and 
characterizing young Jovian planets. SPHERE’S overall scientific capabilities are broadly 
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similar to those of GPI, and it has been operational since mid-2014. A third project, with 
similar goals, and with a high angular resolution that should permit detection of planets at 


r ~ 4AU, is the Subaru Coronagraphic Extreme-AO Project, (Martinache & Guyon 2009), 
which is scheduled to begin scientific operations in 2015. 

Several other methods can also be used to detect extrasolar planets. The precise timing 
of eclipses of eclipsing binary stars has the potential of revealing the masses and orbits of 
unseen companions, and indeed, such timing measurements have been used to conhrm the 
planetary nature of the transiting circumbinary planet Kepler-16 b ( Doyle et al.j2011 ). Spec¬ 
troscopy could be used to identify gases that would be stable in planetary atmospheres but 
not in stars, and Doppler variations of such signals could yield planetary orbital parameters. 


As an example, Brogi et ah (2013) have obtained tentative detections of GO and H2O in the 


emission spectrum of 51 Peg during superior conjunction of its non-transiting short-period 
planetary companion. Starlight scattered from a planetary atmosphere will be polarized (for 


detail, see Seager et ah 2000), and the expected polarization signal lies at the threshold of 


current detectability (Berdyugina et ah 2008 Wiktorowicz 2009). Radio emissions similar 


to those detected from Jupiter could reveal the presence of extrasolar planets, although to 


date, all attempts at detection via this channel have been unsuccessful, see, e.g., Lazio et ah 


(2010). Finally, the detection of artihcial signals from an alien civilization would potentially 


imply the presence of a habitable terrestrial-sized planet or giant planet satellite. 

The Kepler spacecraft has broken through the psychologically important threshold of 
one Earth radius, identifying numerous Earth-size candidates, and has validated several 
planets smaller than Earth, including Kepler-37 b, a planet about the size of Earth’s Moon 


(Barclay et ah 2013). The Doppler radial velocity technique is approaching the comparable 
1 M0 threshold, with the Geneva Planet Search Team’s detection of an Msini ~ 1.1 M0 
planet in a short-period orbit about Alpha Gentauri B (Dumusque et ah 2012). All of 


the Earth-size exoplanets verihed thus far, however, have orbital periods of a few weeks 
or less, and are not worlds of genuine astrobiological interest. Nevertheless, it has been 
unambiguously established that planets of masses somewhat greater than Earth’s mass and 
larger are extremely common in the galaxy, and the emphasis of the held is shifting from the 
detection of extrasolar planets to their physical characterization. In the same manner that 
a new branch of astronomy has been created, a new and important subheld of planetary 
geophysics has also come into being. 
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4. An Overview of the Current Planetary Catalog and Census 


Including both verified exoplanets and high-quality Kepler-detected transiting planet 
candidates, thousands of planets have now been identified, and the outlines of the galactic 
planetary census are emerging. Much of the current inventory of planets is shown in Figure 2. 
In this diagram, which combines planets and planet candidates from several sources, the log 
of the mass ratio of the planet to the central star, log^Q(Mp/M*), is plotted on the ordinate, 
and log^o P is charted on the abscissa. For the case of Doppler detections, we confiate Mp 
with Mpsin(i), and for high-quality candidate planets detected photometrically, we have 
transformed from planetary radius, Rp to Mp using a simple transformation 


Mp= 1 

(u) 

Rp < R-e; 

(2) 

Mp= 1 


Ri > Rp > Re) 

(3) 

Afp — Mj 

Rj > Re- 

(4) 


In the above, the relation ( Lissauer et al.j2011a ) provides a best¬ 

fitting power-law to the Solar System planets spanned in mass by Earth and Saturn. That 
such a simple power law does such a good job in the face of different formation mechanisms 
for the solar system planets can largely be attributed to the fact that only four planets are 
involved in the fit, and to the substantial similarity between Uranus and Neptune. The 
Mp = Mj relation is applied to just a handful of candidates that appear in multiple-planet 
systems, and does not affect the overall trends exhibited by the diagram. In the absence of 
solid information, the relation Mp = (i?p/R 0 )^ M 0 simply assumes that the small number 
of planet candidates with radii smaller than Earth have Earth-like densities. As we discuss 
below, the results from planets that have had both Mp and Rp independently measured 
indicate that there is no single planetary mass-radius relation, and that the dispersion in 
radii for planets of given mass is substantially larger than was generally expected. 

The degree of incompleteness that characterizes Figure 2 is somewhat difficult to eval¬ 
uate, especially in light of the fact that it draws together planets detected by a variety 
of methods. Nevertheless, a reasonable estimate for the current limiting sensitivity of the 
Doppler Surveys, in the absence of highly-focused, high-cadence campaigns on bright nearby 
stars, is A = 2 m s“^, multiplied by an additional period-dependent term of order unity that 
favors shorter orbits. A curious, and mostly unpredicted feature of the diagram is that the 
planets fall into three groups. Observational completeness, furthermore, while not perfectly 


known, see, e.g.. Gumming et ah (2008), is sufficient to ensure that the gaps separating the 


three distinct populations represent real features of the galactic planetary census. 
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Centered on Mp/M* ~ 1 x 10“^ (i.e., ~ 1 Mj), and P ~ 3days, is a population of hot 
Jupiters. The members of this group, on account of their ready detectability, are the best- 
studied (if still incompletely understood) population of extrasolar planets. In aggregate, hot 
Jupiters are clearly gas giants dominated by hydrogen - helium compositions, but in nearly 
every observable property - radii, chemistry, emission properties, density, orbital elements, 
and so forth, they display a tremendous variability that dehes easy categorization. 


Hot Jupiters constitute a prominent fraction of the planets in Figure 1 because they are 
readily detectable, but as a whole, they represent a small fraction of the entire population 
of planets. Among a large sample of main-sequence host stars, they are the exception rather 
than the rule, with an occurrence rate among nearby Sun-like main-sequence stars of only 
about ~1% ( Mayor et al.||2011a ; Wright et al. 2012). From the Batalha et ah (2013) catalog 
of Kepler candidate planets, one estimates a hot Jupiter occurrence rate of 0.5%, in contrast 
to the order-unity occurrence rate of smaller planets with orbital periods P < 100 d. 
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Fig. 2.— Green circles: log;^o(^sateiiite/A^primary) and log]^Q(P) for 634 planets securely 
detected by the radial velocity method, including those hrst identihed via transits. Red 
circles: log]^o(^sateiiite/^primary) and log;^o(-P) regular satellites of the Jovian planets 

in the Solar System, along with Earth’s Moon. Gray circles: Iogio(d4sateiiite/d4^primary) and 
log^o(-P) 1501 Kepler candidates and objects of interest in which multiple transiting 


candidate planets are associated with a single primary. Radii, as reported in Batalha et ah 


(2013), are converted to masses using the mass-radius relations described in the text. Black 


circles: Planets detected by direct imaging, with mass ratios estimated from infrared emission 
detected from the planetary candidate (imaged planets, e.g. HR 8799 b, c, d, & e with 
P > 10"^ d are not shown). Gold circles: Planets detected via microlensing. Purple circles: 
Planets detected via pulsar and white dwarf timing. Blue circles: Jupiter, Earth and Venus. 
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The upper-right hand corner of Figure 2 shows a second population of planets that have 
been revealed by the Doppler velocity surveys. In the absence of a consensus terminology, 
we refer to this gronp as the Eccentric Giants. Planets in this gronp are characterized by (1) 
orbital periods in the range of hnndreds to thousands of days, (2) masses that are, on average, 
several times larger than the mass of Jnpiter, with a median value of Msini = 1.94 Mj, 
although the mean of the nnderlying distribntion is clearly smaller becanse more massive 
planets are easier to detect, and (3) typically snbstantial eccentricities, with a median valne 
e = 0.22, and e > 0.5 not nncommon, although a few members travel on nearly circular 


orbits. The aforementioned HD 114762 b (Latham et ah 1989), with Msini = 11 Mj, and 


P = 114 d is a member of this population. Of order P ~ 5% of the solar-type F, G & K 
main-sequence dwarfs in the local galactic neighborhood harbor an eccentric giant, making 
them roughly ten times more common than hot Jupiters. Like the hot Jupiters, however, 
the occnrrence rate of eccentric giants depends strongly on the metallicity of the parent star 
( Gonzale'z||1997 , 1999). This strong correlation imparts an important, albeit still imperfectly 
nnderstood, cine regarding the natnre of the planet-forming process. The most metal-rich 
stars, with [Fe/H] > 0.3, i.e., more than twice the solar value, are a factor of hve to ten 
times more likely to harbor a hot Jnpiter or an eccentric giant than are stars with half a 
solar metallicity (i.e., those with [Fe/H] < —0.3) ( Santos et al.|2001 ; Fischer fc Valenti||2005 ). 
Stars with M* > Mq are also statistically more likely to harbor a hot Jupiter or an eccentric 
giant. 

Jupiter lies near the long-period edge of the currently known population of eccentric 
giants. Until recently, this placement near the tail of the distribntion was attribnted to 
observational incompleteness among genninely jovian analogs with Jnpiter-like masses and 
Jnpiter-like periods. That is, the apparent dearth of planets with P ~ 10 years and M < Mj 
was thonght to be the simple consequence of long reqnired time scales for observation, and 
the difficulty in detecting low-amplitnde signals. Such assumptions, however, are increas¬ 
ingly being called into question. As an example, consider the stars in the Keck radial velocity 


survey as maintained by the Lick-Garnegie Exoplanet Search (Vogt et ah 2010). By the end 


of 2012, the Keck I telescope had been used to obtain 31,910 Doppler measurements in 
the overall radial velocity database (binned to a minimum cadence of two honrs in the fre- 
qnently occnrring case of mnltiple back-to-back exposures that average over stellar p-modes) 
on 1,298 individnal, primarily solar-type, main seqnence dwarfs. The median derived in- 
strnmental nncertainty on these measurements was a = 1.45 ms“^. Stellar jitter, depending 
on the star, generally contribntes cr = 1 — 5ms“^ in quadratnre to the instrumental nncer¬ 
tainty, and for observations associated with individnal stars, the median time base line was 
HJDobs.iast — HJDobs.first= 2,426 days. Fnrthermore, there were 237 stars with time baselines 
longer than Jupiter’s 4333 d orbital period. The stars in this long-baseline gronp tend to 
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be similar to the Sun in terms of basic properties, and are generally chromospherically quiet 
and photometrically stable, averaging 55 Doppler measurements each. A true-Jupiter ana¬ 
log (P = 11.862 years, K = 12(sini)m s“^) orbiting almost any of these stars would be 
immediately evident. 

Of potentially equal signihcance to the concentrations of planets in the mass - period- 
ratio diagram are the regions of parameter space where the number density of worlds is low. 
There is some evidence that Nature tends to avoid siting giant planets with orbital periods in 
the 10 d < P < 100 d range; the underlying cause of these broad areas of lower concentration 
is not fully understood. There appears to be a general dearth of planets of both short and 
intermediate period with masses ranging from roughly 40 M0 to 100 Mq, consistent with the 
predictions of rapid growth through this mass range in the core-nucleated accretion model 


of giant planet formation (Pollack et ah 1996). 


Perhaps the most important exoplanet-related discovery of the past ten years has been 
the realization that of order half of the Sun-like stars in the solar neighborhood are accom¬ 
panied by systems of one or more planets with periods ranging from days to months, and 
masses falling in the 1 M0 < Mp < 50M0 range (Mayor et ah 2009, 2011b). These plan¬ 


ets are clearly visible as a third population that is distinct from the hot Jupiters and from 
the eccentric giants. Its members constitute the bulk of the transiting planet candidates 
in multiple-transit systems detected by Kepler, which are indicated in gray on the mass - 
period-ratio diagram, and they tend to have relatively low orbital eccentricities. There is 
clearly a desert (a dearth of objects), between giant planets, and these “ungiants” - the term 
of art that we will employ for what is largely still terra incognito - the planetary mass range 
M 0 < Mp < 30 M0. 

Whereas the period distribution of giant planets reveals that hot Jupiters are a separate 
population, distinct from the more distant giants, whose abundance peters out at periods 
of less than a few months, no such distinction is seen among the ungiants. The population 
of these planets appears roughly constant in log period space (as would be expected from 
stability limits imposed by planetary dynamics) for orbital periods between 10 and 100 days. 
Data are too sparse to extend this conclusion to longer periods. At shorter periods, there is 
clearly a drop in the abundance of ungiant planets, but this is gradual, with no counterpart 
to the spike in the 3-4 day range seen for giants. Neptune-size planets are detected down to 
orbital periods of 2 days, and planets with Rp < 2 R0 are observed at periods well under 
one day. The current record-holder for ultra-short period planets is Kepler-78 b, a refractory 
world (it has been detected both in transit and using RV) slightly larger than the Earth with 


an orbital period of only 8.52 hours (Sanchis-Ojeda et ah 2013). 


Little is known about the physical properties of ungiants, apart from that by mass they 













are composed primarily of elements heavier than helium, in contrast to giant planets, and 
that they are a diverse population in terms of dominant constituent by volume. It is not 
clear, for example, whether most of its members are properly described as “super-Earths” or 
“sub-Neptunes” or something else entirely. Curiously, our own Solar System, which contains 
nothing interior to Mercury’s 88-day orbit, lacks any of the larger and closer-in planets that 
dominate this part of the exoplanet diagram; how common terrestrial planet conhgurations 
similar to our Solar System are is unknown, although many stars clearly host different and 
easier-to-detect systems of ungiants. 
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Fig. 3.— Mass-Radius diagram for Extrasolar planets. This figure shows planets for which 
both mass and radius estimates exist, with orbital period coded to run from one day (blue) 
through 100 days (red). Three sample mass-radius relations are shown. The short-dotted line 


is drawn from a simple in-situ formation theory as described in Chiang & Laughlin (2013), 


and is based on rock-iron planets surrounded by hydrogen-helium envelopes, and does not 
include the possible effects of insolation-driven atmospheric mass loss. The long-dashed line 
is the best-£t Mp/M0 = (Rp/Rq)^-^® power-law appropriate to the Solar System planets 
spanning masses from Earth through Saturn. The medium-dashed line is the dMp/M^ = 
Rp/R0 relation derived by Lithwick et al. (2012) from consideration of the transit timing 


variations observed by Kepler for 22 multiply-transiting candidate planet pairs in the vicinity 
of mean motion resonances. 
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A confounding factor is presented by the difficulty in obtaining precise Doppler radial 
velocity-based or transit timing-based mass estimates for super-Earth/sub-Neptune class 
planets. Figure charts a selection of planets for which both mass and radius (and hence 
density) measurements have been made. Even with the log-log scaling, the plot indicates 
to the eye that planets are subject to a wide variety of structures. For example, to within 
the errors, it appears that plants of M ~ 6 M 0 vary in radii by a factor of at least three. 
When mass is viewed as a function of radius, the transition in composition appears sharper, 
with planets i?p > 2 R 0 generally not dominated by rock (including iron), and those with 
radii Rp <1.6 R 0 and measured masses mostly rocky. While the upper limit to the size of 
abundant rocky planets within a ~ 0.5 AU of their stars seems well-established, it is quite 
possible that many or even most planets smaller than 1.6 R 0 are composed primarily of 
low density constituents, but that the low masses of such planets makes them difficult to 
measure. 


On Figure we have over-plotted three potential mass-radius relations for ungiants. 
The first relation is the Mp/M^ = best-fit to the Solar System planets spanned 

in mass by Earth and Saturn. Despite the uncertainties, and the clear intrinsic variations, 
this relation provides a reasonable £t to the mass-radius relation of planets smaller in mass 
than Saturn. The Solar System £t is bracketed by two alternate relations. The first, with 
a steeper slope, is derived from the assumption that low-mass, short-period planets have 
rocky cores surrounded by hydrogen-helium envelopes, with envelope masses dictated by in 
situ formation, as discussed in Chiang & Laughlin (2013). The second alternate relation, 
3 Mp/M 0 = i?p/R 0 , has a shallower slope, and is a fit derived from an ensemble of transit 
timing measurements applied to 22 systems within Kepler’s multi-transit transiting planet 
candidate population by Lithwick et al. (2012). Recent work by Weiss & Marcy (2014) 
has generated additional Doppler velocity-based mass measurements for a selection of Ke¬ 
pler candidates that are not shown in Figure]^ These latest results reinforce the conclusion 
that planets with Rp < 4 R 0 span a very wide range of masses. 


Kepler is capable of detecting rocky terrestrial planets analogous to Venus and Earth, 
but the spacecraft is unlikely to hud a statistically robust sample of terrestrial planets similar 
to those in our Solar System. The mass-period diagram is thus still largely incomplete for 
terrestrial planets with Solar System-like masses and orbits, and will likely remain so for 
some time. In the absence of solid data, we can only resort (with an abundance of cautionary 
disclaimers) to statistical extrapolations to anticipate the inventory of low-mass planets. 


To construct a concrete (albeit uncertain) extrapolation, we assume that planet for¬ 
mation is highly efhcient, and that long-distance migration (Aa > a), where a is the cur¬ 
rent semi-major axis of the planet and Aa the amount it has migrated, in the super-Earth 
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population does not regularly occur. We then idealize the population of Kepler planetary 
candidates from Batalha et ah (2013) as being composed of solar-composition solids, and we 
assign each planet a surface density, cxsoiid = Mp/27rap, where Op is the orbital semi-major 
axis, computed with the approximation that M* = Mq, and Mp = M®. The 

resulting 7V=1,925 surface densities obtained for low-mass candidates with R < 5 Rq are 
plotted in Figure]^ The best-fit power law to the median data can be interpreted to define a 


“minimum mass extrasolar nebula” Usoiid = 50 (a/1 AU)~^'^gcm“^, which has a very similar 
power-law index, a = —1.6 to the canonical minimum mass solar nebula (MMSN) value of 
a = —1.5 (Weidenschilling 1977), but with a density normalization that is 5x larger. 



Fig. 4.— The solid surface density profile of the “Minimum-Mass Extrasolar Nebula” 
(MMEN) constructed from Kepler data Chiang & Laughlin (2013). Gray circles: cTsoiid.i = 
Mi/2na‘f, computed from Rep/er planets with Rp < 5 Re and P < 100 d, assuming Mp/Me = 
(Rp/Re)^'*^® and solar-mass host stars. For reference, Usoiid for the minimum-mass Solar 
nebula (MMSN) — or technically its extrapolation inward, since no planet is present interior 
to Mercury at a = 0.4 AU — is plotted as a blue solid curve. As plotted, the MMEN is 
a factor of 5 times more massive than the MMSN at these distances. Qgas < 1, where Q 
is the Toomre stability criterion (|Toomre 1964), implies the nebular surface density in gas 


required for axisymmetric gravitational instability. 


Extending the Usoiid = 50 (a/1 AU) ^'^gcm ^ power law to a distance a ~ 1 AU leads 
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to two suggestions, or better, to two moderately informed speculations: 

• Roughly half of the sunlike stars in the Sun’s neighborhood will harbor a planet as 
large as (or larger than) Earth in (or near) an Earth-like orbit. 

• Among the planetary systems of single main-sequence stars, the inventory of mass in 
the Solar System’s terrestrial planets will lie below the median total planetary mass 
interior to a 2 AU. 


5. Planetary Geophysics 


Investigations of extrasolar planets can be divided conceptually into (1) analyses of 
planetary structure and composition, (2) the exploration of planetary orbital dynamics, 
and (3) the study of planet formation. These three areas closely inform one another. For 
example, at a given planetary age, Tp, the radius of a giant planet (a zeroth-order structural 
property) depends on the mode of formation. A planet that emerges from the process of core 
accretion will probably be smaller at early Tp than one that arises from disk gravitational 


instability (Marley et al. 2007). The planetary radius, in turn, has a strong influence on the 


tidal dissipation rate, with tidal luminosity obeying E oc (i?p/a)^, and frequently, the long¬ 
term dynamical behavior of a compact planetary system is sensitive to its tidal evolution, 
e.g., Mardling & Lin (2002). In the remainder of this article, we will give an overview of 
the “geophysics” of extrasolar planets, with an eye toward elucidating the connections to 
theories of dynamical evolution and theories of planetary formation. 

The direct characterization of extrasolar planets has been enabled in large part by elab¬ 
orations on transit photometry, and as a result, our understanding is strongly biased toward 
short-period planets orbiting close to their parent stars. Most of the well-characterized gi¬ 
ant exoplanets were discovered by ground-based transit searches, with most of the ungiants 


found by the Kepler spacecraft (Lissauer et al. 2014a). With this context, we then provide 


a synopsis of the theory of planet formation, largely informed by the contents of our own 
Solar System. 


5.1. Exoplanetary Systems 

Radial velocity surveys have detected scores of systems with more than one exoplanet, 
including more than one dozen nearby stars with three or more planets. In systems with high 
multiplicity, the planets tend to have low masse and low eccentricities, although a number 
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of systems with three or more giant planets are known. Planets tend to be clnstered into 
systems in the sense that there are more mnltiple planet systems than wonld be the case 
were planets randomly distribnted among stars. Hot Jnpiters are an exception to this rnle 
- few of them are known to have companion planets. (Nonetheless, the star v Andromedae 
has at least three jovian-mass planets, one of which is a hot Jnpiter. The remaining planets 
are far more distant from v And and travel on eccentric orbits. 


The mnlti-transiting systems from Kepler provide a mnch larger and very rich dataset 
of short-period (primarily from a few days to several months) planets that can be nsed to 
powerfnlly test theoretical predictions of the formation and evolntion of planetary systems. 


The major hndings, npdated from (Lissaner et ah 2011b) are: 


1. The large nnmber of candidate mnltiple transiting planet systems observed by Ke¬ 
pler show that flat (i.e., nearly coplanar) mnlti-planet systems are common in short- 
period orbits aronnd other stars. This resnlt holds for nngiant planets in the size range 
of i?p ~ 1 - 6 Re, bnt not for giant planets. Not enongh data are yet available to 
assess its viability for smaller worlds, nor for planets with orbital periods longer than 
a few months. 

2. Mean-motion orbital resonances and near resonances are mildly statistically favored 
over random period ratios, bnt most planets are neither in nor very near snch res¬ 
onances. First-order resonances (e.g., period ratios of 2:1, 3:2, etc.) dominate, bnt 
second-order resonances (with period ratios 3:1, 5:3, etc.) also are manifested. Most 
near-resonant planet pairs have period ratios from one to a few percent larger than 
those of the nearby resonance. Few planet pairs have period ratios slightly smaller 
than those of hrst-order resonances. 

3. Almost all candidate systems snrvive long-term dynamical integrations that assnme 
circnlar, planar orbits and the mass-radins relationship given in eqnations (2) - (4), 
derived from planets within onr Solar System. 

4. Attempts to generate a simnlated ensemble of planetary systems to match the ob¬ 
served Kepler planets do not strongly prefer hts by a single homogeneons popnlation. 
In particnlar, there is some evidence for a popnlation or snbpopnlation of systems that 
either contains only one detectable planet per star or mnltiple planets with high rela¬ 
tive orbital inclinations that prodnces the observed excess of singly-transiting systems. 
A single popnlation, whose characteristics are described below, can acconnt for the 
vast majority of mnlti-transiting systems, as well as ronghly half of the single planets 
observed. A third, rarer, gronp of planar densely-packed mnlti-planet systems, also 
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appears to be present. Note that these populations of planetary systems need not be 
cleanly separated, i.e., there may also be signihcant numbers of intermediate systems. 

5. The abundance and distribution of multiply transiting systems implies that the mean 
number of planets in the 1.5 R® < i?p < 6 R® and 3 < P < 125 day range per 
star for those stars with at least two such planets of this radius is ~3. According to 
Kepler observations, only ~5% of stars have planetary systems within these size-period 
limits. 


6 . The inclination dispersion of most candidate systems with two or more transiting 
planets appears to have a median of ~ 2° (e.g., Lissauer et ah 2011b), suggesting 
relatively low mutual inclinations similar to the Solar System. 


7. Many Kepler targets with 1 or 2 identihed planet candidates must be multi-planet 
systems where additional planets are present that are not transiting and/or too small 
to be detected. 


5.2. The Radii, Densities, and Compositional Structures of Extrasolar Planets 

Hot jupiters, which have P < 7 d and Mp > 100 Mq, represent extrasolar planets at 
their most alien, but they are nevertheless by far the best studied. The short orbital periods 
and large masses of these worlds allow for precise and full orbital characterization for cases in 
which transits occurand their high surface temperatures generate infrared planetary fluxes 
that are readily detectable. Scores of short-period giant planets now have accurate physical 
measurements, but the spectra and even the radii of these planets have dehed straightforward 
explanation. The origins, furthermore, of the majority of the physically characterizable 
extrasolar planets are not yet fully clear. It is generally thought that conditions at a < 
0.1 — lAU in protostellar disks are not favorable to the successful completion of the core 


accretion process ( 

Bodenheimer et ah 2000 

). It is therefore believed that planetary migration. 

either via quiescent disk processes ( 

Lin et al.f 

1996; 

Lubow & Ida 

2011 

), via planet-planet 


scattering (Juric fc Tremaine||2008 ) or via Lidov-Kozai Cycling with Tidal Friction (Eggleton 


& Kiseleva-Eggleton 2001 Fabrycky & Tremaine 2007), or some combination thereof, plays 


an important role in delivering the planets to their current locations. 


®Transit probabilities are given by Vtr = 0.0045 


/ \ 

1+e cos(7r/2—cj) 

V R© ) 

l-e2 


, where w is the longitude 

of periastron of the star such that uj = 90° corresponds to the prospective transit midpoint, and e is the 
orbital eccentricity. V ~ 10% for typical hot Jupiters. For a distant observer sited at a random viewing 
angle, Re = 0.5%. 
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Fig. 5.— Densities and masses for planets that have individnally measured values for and 
Mp. Theoretical mass-radius relations for planets of various pure compositions are shown for 
comparison. Note that many of the measurements for low-mass planets were derived from 
transit timing variations, and these data points tend to have large uncertainties. The six 
most massive Solar System planets (Venus, Earth, Uranus, Neptune, Saturn, and Jupiter) 
are included on the diagram, and can be distinguished on account of their low effective 
temperatures. 


The characterization of extrasolar planets with Mp < 30 M 0 is difficult, even when 
the orbital periods are short. Low mass planets have small transit depths, small radial 
velocity half-amplitudes, and generally small transit timing variations, which make it hard to 
accurately determine densities. Nevertheless, dedicated follow-up on small-radius transiting 


planets identihed by the Kepler Mission has begun to yield results (Weiss & Marcy 2014 


Jontof-Hutter et ah 2014), and Figure shows that the broad outlines of the planetary 
distribution for M® < Mp < 30 M 0 are being delineated. While this diagram is generally 
consistent with the p{Rp) values exhibited by the Solar System planets, it is clear that the 
ungiant regime spanning (M 0 < Mp < 30 M 0 ) encompasses a very wide range of planetary 
densities (and by extension, planetary compositions). Frustratingly, there is substantial 
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compositional degeneracy if planetary mass and radius are the only quantities that are 
measured. Often, a substantial range of admixtures of metallic, rocky, icy, and enveloping 
gaseous components can jointly contribute to produce a planet of given mass and radius 


(Adams et ah 2008). The lack of significant concrete physical data for low-mass planets 


means that guidance must be drawn from theoretical work. The internal structures of planets 


more massive than Earth are reviewed in Baraffe et ah (2014) and the references therein. 


The expected properties of “ocean” planets with large water fractions have been discussed 
by a number of authors, see, e.g., (Sotin et ah 2007| 2011). 

Figure suggests that the overall distribution of planetary densities reaches a broad 
minimum in the region pmin ~ 0.2 —>■ l.Ogcm”^, at roughly Mp ~ 0.1 Mj. Above this mass, 
planetary compositions are dominated by hydrogen and helium. Over 130 well-characterized 
examples of transiting H-He dominated giant planets (with Mp/Mq > 100) are now known. 
It is clear that within this higher mass regime, variations in mass and stellar insolation are not 
the only factors responsible for the observed range in planetary sizes, and radius anomalies 
- observed planetary sizes that are substantially different from those predicted by structural 
models - have been evident since the discovery of the first transiting extrasolar planets. 


Early structural evolutionary models for giant planets computed by, e.g., Bodenheimer et ah 


(2001) and Guillot & Showman (2002), suggested that a H-He dominated planet with HD 
209458b’s mass, insolation, and age should have a radius of approximately ~ 1-lRj, a 


hgure that is startlingly at odds with the observed value, = 1.38 ± 0.02 Rj (Southworth 


2010 ). 


During the last decade, the accumulating measurements of giant planets have demon¬ 


strated that radius anomalies are by no means anomalous, see, e.g.. Burrows & Orton (2011), 


or 


Baraffe et al. (2010) for overviews. The substantial variation in observed radii for the giant 


planet population is indicated by Figure]^ which charts the measured radii and uncertainties 
for the transiting giant planets with accurately measured masses against an insolation-based 
orbit-averaged planetary equilibrium temperature given by 


1/2 


T* 


\2a) ( 1 - 02 ) 1 / 8 ’ 


(5) 


which implicitly assumes that the planet has an albedo of zero and re-radiates energy uni¬ 
formly over its entire surface and depends on the average flux received by the planet over 
its orbit. Incorporating the effect of albedo results in multiplying the right-hand side by 
(1 — A)^/^. In Equation]^ T* is the stellar effective temperature, e is the planetary orbital 
eccentricity, i?* is the parent star’s radius, and a is the planetary semi-major axis. The ac¬ 
tual temperatures in most of the exoplanetary atmospheres are unmeasured, and, depending 
on the atmospheric properties, chemistry and dynamics, may be higher or lower than Tgg. 
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Fig. 6.— Planetary radii, i?p/Rj, versus insolation-derived planetary equilibrium temper¬ 
atures, Teq (see Equation]^, for 130 planets with Mp/M0 > 100 that have both transit 
and radial velocity characterizations. The planets are color-coded by mass, suggesting that 
planets substantially more massive than Jupiter are less likely to exhibit signihcantly inflated 
radii. 


At a given planetary mass, it is expected that the radius of a mature gas-giant planet 
should primarily be determined by the amount of radiative energy that it receives from its 
star. This expectation is supported by the clear trend in Figure]^ which indicates that Rp 
is positively correlated with planetary Teq. The degree of quantitative agreement can be 
obtained by computing model radii, Rm, for H-He composition planets as a function of Mp 
and Teq, and comparing with those of the observed planets. We dehne the radius anomaly, 
TZ as the difference 71 = Rohs — Rpved between the observed and predicted radii. Figure 
shows the result of such a comparison to a set of baseline core-free models computed by 
Bodenheimer et ah] (|2003) using a Henyey-type planetary structure and evolution code. 


Figure [^indicates that the structural models do explain some, but not all, of the observed 
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Fig. 7.— “Radius anomalies”, 77 = Robs — Rpred, versus insolation-based planetary equi¬ 
librium temperature, Teq. Planetary models which account for insolation typically predict 
radii that are too small when Tgq > 1200 K, and too large when Tgq < 1000 K. 


variation in planetary radii. For Tgq > 1000 K, planets are usually larger (and sometimes 
substantially larger) than predicted by the models. At lower temperatures, radii tend to be 
smaller than what one would expect for a giant planet of solar composition. Several authors, 
Enoch et ah (2011), have noted the trend that is evident in Figure]^ and a power-law 


e.g. 


fit to the current data suggests an 77 oc dependence. 

Giant planets with negative 77 probably contain a substantial fraction of heavy elements. 
The majority of observed planets with Tgq < 1000 K have negative radius anomalies; for this 
population TZ = —0.12, which suggests that the average short-period giant planet has a 
Mcore ~ 30 Mq core mass. Planets orbiting metal-rich stars are statistically associated with 
larger core masses (Miller & Fortney 2011), and there is likely an inverse relation between 
planetary mass, Mp, and the degree of a planet’s enrichment in heavy elements, Zp/Z^ 
(Laughlin et al. 2011). 
































































In order to inflate an H-He dominated planet to large TZ, models must appeal to a 
cryptic source of heating at adiabatic depth, and a wide range of heating mechanisms have 
been proposed. Those that have received significant attention include dissipation stemming 


from (1) tidal orbital circularization (Bodenheimer et al. 2001), (2) obliquity tides (Winn 


et al. 2005), but see the cautions raised by Levrard et al. (2007), and see Fabrycky et al. 


(2007) for full-blown arguments against obliquity tides, (3) “kinetic heating” in which wind 
energy is converted into heat ( |Guillot fc Showman||2002 ), (4) dissipation induced by thermal 
tides (Arras & Socrates 2010), and (5) Ohmic dissipation (Batygin & Stevenson 2010). 


A clue is provided by the TZ oc dependence. It is difficult to understand this 
correlation if the radius anomalies stem primarily from tidal heating, for which one would 
expect the major dependence to be with eccentricity, 7Z = TZ{e). Individual planets, of 
course, may be significantly affected by dissipation, and there is a wide observed dispersion 
in planetary size at all Teq. 


Arras & Socrates (2010) revived the thermal tidal torque mechanism that was introduced 


by Gold & Soter (1969) to account for Venus’ low-frequency retrograde spin. In Arras & 


Socrates (2010 )’s picture, uneven global heating of a short-period planet produces a compar¬ 


atively cold, high-density bulge that trails the star-planet line. Gravitational torque acting 
on the bulge spins the planet up into a state of asynchronous rotation, enabling steady-state 


constant-lag tidal energy dissipation to occur in the planet. Arras & Socrates (2010) suggest 


that this dissipation may be sufficient to inflate hot Jupiters to their observed radii. Other 


workers, however, e.g., Gu & Ogilvie (2009), argue that rapid isostatic re-adjustment of the 


planet will offset the thermal bulges for forcing frequencies (e.g., the planetary mean mo¬ 
tion, n = 271/P) that are lower than the planetary dynamical frequency, which is of order 
the free-fall time scale, 


The kinetic heating mechanism, as outlined by Guillot & Showman (2002), posits that 
a constant fraction, rj ~ 0.01, of the total flux received by a planet is converted into kinetic 
energy by atmospheric gradients, and that this energy is dissipated at adiabatic depth within 
the planet. While the flow patterns on the surfaces of strongly irradiated extrasolar planets 
are both unknown and a matter of debate (see below), an overall E oc scaling is implied 
by the kinetic heating hypothesis, leading to 77. oc which is somewhat weaker than the 
observed 77 oc dependence. As was the case with the evaluation of the viability of tidal 
heating, this argument does not imply that kinetic heating is absent, but rather, that it is 
unlikely the primary cause of the radius anomalies. 


Batygin & Stevenson ( 2010| argue that ohmic dissipation inflates hot Jupiters. In 
their model, significant atmospheric electrical conductivity at Teq > 1200 K is enabled by 
valence electron ionization of alkali metals. The presence of Na in the atmosphere of HD 
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209458b has been directly observed via transit spectroscopy ( Charbonnean et al.||2002 ). The 
Batygin-Stevenson model assnmes that the atmospheric layers of the swollen planet have 
rotation-indnced zonal jets, which are seen on all fonr Solar System giants, and which are 
a generic featnre in two-and three-dimensional hydrodynamic simulations of exoplanetary 


surface flows (Showman & Polvani 2011). Conductive atmospheric zonal flows introduce a 


meridional surface current that, in the Batygin-Stevenson model, flows from the planetary 
poles to the equator, and then closes the circuit by returning poleward through the deep 
convective interior of the planet. Resistive heating within the planet, especially in the near¬ 


surface (but still adiabatic) layers, provides an internal energy source. Batygin et ah (2011) 


hnd that this heating is sufficient, for reasonable planetary atmospheric compositions and 
magnetic held values, to provide 71 values of order those that are observed. 

If one further assumes that the planetary Elsasser number, A = aiB'^/2pQ, is roughly 
constant across the aggregate of transiting planets, the Batygin-Stevenson theory implies an 
71 oc dependence. This power-law index a = 2.7 is steeper than the observed value 
of a ~ 1.4. A more sophisticated treatment, however, must account for the back reaction 
onto the wind velocity, v, generated by the Lorentz force. Three-dimensional simulations 


by Rauscher & Menou (2013) suggest that the Lorentz force provides a drag term to the 


velocity held, which in turn reduces the current and decreases the dissipation rate, driving 
the power-law index a to lower values. In summary, it appears likely that Ohmic dissipation 
may provide the bulk of the explanation for the large radius anomalies that characterize 
many short-period hot Jupiter type planets. 

The realm of the gas giants has a lower mass bound of about 50 M®. Figure (which 
charts masses versus densities) suggests that for Mp < 50 M 0 , the overall fraction of refrac¬ 
tory material increases steadily toward lower masses. This observed trend is in part due to 
selection ehects - the lowest mass planets for which both R and M can be measured have 
short periods and therefore, very high surface temperatures. As an example, Kepler-78 b 
( Howard et al.||2013 Pepe et ah||2013 ) has Mp ~ 1.7 M 0 , Rp ~ 1.1 R 0 , and pp ~ 5.5gcm“^, 
and must therefore be composed largely of metals and silicates, but its orbital period is only 
P = 0.35 days, and it is so close to its star that it would very likely have lost any low-density 
component to its original mass. Indeed, the low-Teq, low-Mp planets shown in Figure 
whose masses have been determined through analysis of transit timing variations, and which 
are thus less sensitive to a bias toward very short periods, have systematically lower densi¬ 
ties than those which have had their masses determined via Doppler velocity measurement 


(Jontof-Hutter et ah 2013). 


Bulk densities give an important clue to the overall planetary structures, but, unfortu¬ 
nately, there are a wide range of compositional structures involving admixtures of metals. 
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rock, ices, and gas that can reprodnce a given planetary mass and radins. See, e.g., Fort¬ 


ney et al. (2007) for a catalog of planetary strnctnral models spanning a broad range of 


compositional assnmptions. 

The tidal Love nnmber, k 2 provides an effective parameterization of the degree of central 
condensation exhibited by a particnlar density strnctnre, p(r), with lower valnes for k 2 


mapping to higher degrees of central condensation. Ragozzine & Wolf (2009) introdnced 
and ontlined the practical problem of measnring fc 2 ’s for individnal extrasolar planets. For 
hot Jnpiters with the shortest periods, they show that the planetary interior indnces apsidal 
precession rates of np to a few degrees per year. 

If a short-period planet in a co-planar, hierarchical (oc 3> Ob), two-planet system has 


tidally evolved to an eccentricity hxed point (see, e.g., Wn & Goldreich 2002 and Mardling 


2007), then the ontcome of this dissipative, deterministic process is that the apsidal lines of 


the two planets are forced to process at the same rate. To a good approximation, when a 
hxed point conhgnration arises, the secnlar gravitational pertnrbation from the inner planet, 
b, fnlly dictates the apsidal precession rate, the, of the onter planet, c. This precession is 
likely too snbtle to be directly observable in a radial velocity time series alone, bnt it is 
nonetheless a fnlly determined fnnetion of directly measnrable orbital parameters, a^, Oc, Cb, 
Cc, Whi Wci and nc, in combination with the masses Mb and M*. The precession rate, thb 


of the inner planet, on the other hand, is dictated by a nnmber of contribntions (Batygin 


et al. 2009), inclnding secnlar pertnrbations from the onter planet, general relativity, and 


the rotational and tidal bnlges of both the inner planet and the star. In the common case 
where the stellar bnlges are negligible, and for an inner planet whose orbital hgnre, mass, 
and radins have all been determined via high-qnality photometric transit and Doppler-based 
observations, ^2 constitntes the only signiheant nnknown. 

As a conseqnence, for systems with favorable orbital conhgnrations, one potentially has 
access to a mechanism that directly probes the interior strnctnre, and which can poten¬ 
tially discriminating between “snper-Earth” and “snb-Neptnne” conhgnrations. To date, 
this method has only been applied to the HAT-P-13 b-c system ( Batygin et al.||2009 ), where 
the inner planet mass is Mb = 0.85 ± 0.038 Mj, well above the Mp < 50 M 0 scale of interest 


here. Kramm et al. (2012) hnd 0.265 < ^2 < 0.379 for HAT-P-13b, which reqnires that any 

icore < 27 Mq. The TESS mission ( |Ricker et al.||2010 ). 


central core in the planet have mass M( 
and CHEOPS missionj^ which are schednled for lannch in 2017, and the PLATO Mission 
(Barban et al. 2013| ), selected in Febrnary 2014 and planned for lannch in 2024, may detect 
nearby transiting mnltiple-planet systems for which the orbital properties can be determined 


'http://Cheops.unibe.ch/ 
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with enough precision to enable accurate measurement of k 2 - Finally, implicit evidence for 
the widespread existence of dense planetary cores is provided by objects such as CoRoT 7b 
(Leger et ah 2009) and Kepler-lOb ( Batalha et ah||2011 ). For a detailed overview of the types 
of systems for which interior structures can be probed, see Becker & Batygin (|2013). 


5.3. Observations of Exoplanetary Atmospheres and their Interpretation 


The foregoing discussion makes it clear that the radii and interior structures of the 
hot Jupiters continue to provide signihcant challenges to a comprehensive theoretical under¬ 
standing. An analogous situation also holds for the atmospheric characterization of these 
planets. 


Nevertheless, despite the extraordinary difficulties inherent in separating the signal of 
a faint planet from a bright star, there have been tangible successes, the majority of which 
have come from space-based measurements. The Spitzer Space Telescope (Werner et ah 


2004) was designed prior to the detection of transiting extrasolar planets, and exoplanetary 


observations were not part of the mission’s scientific requirements. Serendipitously, however, 
Spitzer has played a key role in probing exoplanetary atmospheres. HST, MOST, CoRoT 
and Kepler have all also made substantial contributions, as have ground-based telescopes 
including Keck and the VLT. The literature on the characterization of the atmospheres 
of short-period planets is now extensive, and space considerations preclude an exhaustive 


summary here. See Seager & Deming (2010), Perryman (2011), Line et ah (2012), Tinetti 


et al. (2013), and Madhusudhan et ah (2014) for successively recent overviews. It is important 


to stress that while the accumulated catalog of measurements is by now impressive, the 
signal-to-noise of individual measurements is often low, systematic errors may, in general, be 
significant, and the number of unknowns, including composition, meteorology and chemistry 
that characterize an atmosphere is daunting. 


Atmospheric characterization has mostly been carried out for planets that transit. Ex¬ 
tant observations divide into (1) occultation studies, (2) transit spectroscopy of varying 
resolution, and (3) elaborations on full-phase spectrophotometry. 

A wavelength-dependent decrease in flux can be measured during occultation (also 
known as secondary eclipse), when a planet passes behind a star as seen from Earthj^ At 
visible wavelengths, the occulted light from the planet consists of reflected and scattered 


®Some planets, such as HD 17156 b, with its highly eccentric e = 0.68, w = 112° orbit, undergo transit, 
but not occultation. 
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starlight, and for hotter atmospheres, emission from the planet itself. The primary con¬ 
clusion is short-period, hot Jupiters are dark. For example, Rowe et al. (2008) used the 
MOST telescope to hnd that the geometric albedo, Ag for HD 209458 is Ag = 0.038 ± 0.045 


in the A = 400 —?■ 700 nm optical bandpass. More recently, Evans et al. (2013) used the 


STIS instrument on HST to obtain an optical and near-UV spectrum (essentially a reflection 
spectrum) for HD 189733b. They measured a geometric albedo, Ag = 0.40 ± 0.12 in the 
A = 290 —?• 450 nm region, and Ag < 0.12 for A = 450 —?• 570 nm (consistent with the Rowe 
et al. result for HD 208458b). The planet’s apparent upturn in reflectivity at short wave¬ 
lengths is suggestive of a cerulean visual appearance, and is attributed to starlight reflecting 
from clouds coupled with high-altitude sodium absorption of red light within the overlying 
clear air. 

At near- and mid-infrared wavelengths, in the Rayleigh-Jeans tail of the star’s emission 
spectrum, a typical hot Jupiter has a flux ratio, Rp/E* > 10“^, and Spitzer has proved to 
be very well suited to detecting exoplanets in occultation. During the cryogenic phase of its 
mission, Spitzer’s IRAC, IRS, and MIPS instruments were used to observe 15 exoplanets in 


occultation, starting with MIPS detection of HD 209458 b at 24 pm (Deming et al. 2005), 


and the IRAC detection of TrES-1 at 4.5 pm and 8 pm ( Charbonneau et al.||2005 ). Following 
depletion of its liquid helium, Spitzer has continued to observe occultations with IRAC at 
3.6 pm and 4.5 pm, and over thirty planets now have eclipse depths measured in both of 
these bands. 

Photometric occultation measurements at a range of different bandpasses amount to a 
low-resolution dayside emission spectrum of an extrasolar planet, and a handful of planets 
have been measured in Eve or more bands, including J (1.2 pm), H (1.6 pm), and K (2.2 pm) 


from the ground, e.g., Croll et al. (2011b), in addition to the Spitzer channels. For the most 


favorably detected transiting planets, such as HD 189733b, Spitzer’s IRS instrument was 
used to obtain mid-infrared dayside emission spectra that were resolved into several dozen 
spectral resolution elements in the 5 pm to 14 pm range, with each spectral element having 
S/N ~ 10, (e.g., Grihmair et al.[ 2008). Also for HD 189733b, HST’s NICMOS has been used 
to obtain near-IR spectra presenting 18 spectral channels in the 1.5 pm to 2.5 pm range, with 


S/N in each element also of order 10 (Swain et al. 2009a). These resolved spectra from IRS 
and NICMOS display signihcant variations that are strongly suggestive of molecular lines 
(including H 2 O, NH4, and CO), yet htted multi-parameter models are imperfect, suggesting 
a substantially incomplete picture of the radiative structure, chemistry, and dynamics of the 
atmosphere. 

Naturally, there has been a vigorous effort to interpret the panoply of occultation results 


with theoretical models. For example, Fortney et al. (2008) suggest that strongly irradiated 
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atmospheres (corresponding ronghly to the half of the observed planets that receive the 
largest orbit-averaged flnxes) have thermal inversions cansed by hardy molecnles snch as 
TiO or VO absorbing and re-radiating starlight at low pressnres high in the atmosphere. 


Knutson et al. (2010) present evidence for an empirical correlation showing that chromo¬ 


spheric activity inhibits the formation of such thermal inversions, perhaps by destroying the 


inversion-producing molecules through an elevated flux of high-energy photons. (Madhusud- 


han 2012), (see also Madhusudhan & Seager 2010), have suggested that C/0 ratios form an 
additional dimension by which planets can exhibit planet-to-planet variation in occultation 
depths. 


Line et al. (2013) discuss a Bayesian retrieval approach - see, e.g. Rodgers (2000) - to 


spectral modeling which uses nonlinear optimal estimation to determine atmospheric tem¬ 
peratures, compositions and structures. This technique quantihes the relative information 
content of a given data set, and provides an estimate of the number of parameters (e.g. dis¬ 
tinct molecular constituents) that the data can usefully be employed to determine. Spectral 
retrieval has been used to address a number of controversies of interpretation, for instance, in 


connection with the apparent presence of methane in the atmosphere of HD 189733b (Swain 


Figure consolidates a full range of broad-band occultation results onto a single plot. 
The planets are ordered according to specihc orbit-averaged received flux, as measured by 
Tgq. Occultation depths in the various infrared bands are transformed to flux ratios, F/ 
relative to the expectation for a zero-albedo dvr blackbody re-radiator. It is clear from the 
diagram that there is little or no systematic degree of similarity among the planets. The 
data are consistent with random flux ratios relative to < iV >= 1.5, with standard devia¬ 
tion, (Ta = 0.5. This distribution of measurements suggests a complex mix of contributing 
factors (stemming either from conditions on the planets, or from errors introduced in the 
measurement pipelines), with the attendant applicability of the central limit theorem. The 
< Fa >= 1-5 mean value, however, implies that on the whole, exoplanets are partially re¬ 
distributing their received flux via re-radiation from the night side. Indeed, Figure is in 


general concordance with the view expressed by Hansen et al. (2014) that, to the precision 


of the measurements, the observed emission spectra of exoplanets are generally consistent 
with blackbodies. 

During transit, starlight flows through the atmosphere of the planet all along the day- 
night terminator. A stellar spectrum obtained during transit thus contains an admixture of 
the planet’s transmission spectrum, which holds out the possibility of obtaining a detailed 
probe of the atmospheric column. An early application of transit spectroscopy was presented 


by Charbonneau et al. (2002), who used the STIS spectrograph on HST to detect sodium in 
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the atmosphere of HD 209458b by noting that the transit depth associated with a band pass 
covering the sodium resonance doublet at 589.3 nm was deeper than in adjacent bands. HST 
has also been used to detect evidence of water vapor, methane, and carbon monoxide (e.g. 
Swain et ah 2009a[b ) in the 1.5 —?• 2.5 /rm transit spectra of HD 209458b and HD 189733b. 


For planets transiting H > 8 primaries, a variety of broad-band transmission spectra 
have been obtained for a variety of planets. Interpretation of these spectra is challenging, and, 
in parallel with the situation for occultations, hrm conclusions are hampered by intrinsically 
low spectral resolution, low signal-to-noise, and the need to adopt a whole range of model 
assumptions. For more detail, the reader is advised to consult the review by [Tinetti et ah 


(2013). JWST is expected to have signihcant utility for transit spectroscopy, and optimism 
exists that it will be able to usefully prove the atmospheres of potentially habitable planets 
transiting M-dwarf stars Deming et ah (2009). 


I » KP ■ J « H « K O 3.6 O 4.5 ■ 5.8 ■ 8.0 [ 



Fig. 8.— A collection of broad-band occultation measurements. Planets are ordered on 
the a:-axis by the equilibrium temperatures, Teq = (i?y^T*)/((2a)^/^(l — that they 

would have if they were Ag = 0 blackbodies re-radiating from the full planetary surface area. 
Measured occultation depths are expressed as ratios of the observed eclipse depth in the 
band of interest relative to the expected depth in that band if the planet were a uniformly 
re-radiating Ag = 0 blackbody. The color coding for each band ass are given in the legend 
and range from Kepler observations in the optical to 8 pm. 
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There is an exciting possibility of directly detecting planetary snrface winds by observ¬ 
ing the Doppler shift of spectral absorption lines that are creating by the starlight shining 


through moving parcels of air at the planet’s day-night terminator. Snellen et ah (2010) used 


the VLT to obtain a high-resolution transit spectroscopy of HD 209458b with R ~ 100, 000 
at wavelengths near A = 2/xm. By cross-correlating on dozens of expected CO lines, these 
authors were able to produce a time-dependent measurement of the average velocity of the 
atmosphere at the planetary terminator during the course of the transit. With the orbital 
motion subtracted out, Snellen et al. (2010) measure a net blueshift of 2 ± 1km s“^, and 


suggest that this represents a global flow from the dayside hemisphere to the nightside hemi¬ 
sphere at atmospheric levels in the range 0.01 bar < P < 0.1 bar. Showman et al. (2013) 


hnd that generic flow patterns with this large net blueshift can be created within current 
extrasolar planet atmospheric circulation simulations (see, e.g.. Showman et al.[ 2009), by 
varying the amount of frictional drag in the atmosphere. 

For planets in the ungiant regime, atmospheric measurements obtained via transit spec¬ 
troscopy can potentially sort out the compositional degeneracies that prevent strong conclu¬ 
sions being drawn from bulk planetary density measurements (e.g., those shown in Figure]^. 
A planet with a refractory core covered by a deep solar-composition atmosphere will have 
a wavelength-dependent transit depth that can be potentially distinguished, using current 
observational sensitivity, from the signature produced by a planet with an atmosphere (and 
an overall structure) dominated by high molecular weight species such as water. For exam¬ 
ple, if the atmosphere of a planet such as Gliese 1214b, with Mp = 6.6 Mq, ppi = 0.35 p0, 
and Tpi 500 K ( [Charbonneau et al.||2009a[ ), is primarily composed of H2, then the resulting 
atmospheric scale height, H = ksT/fig, will be of order H ~ 200 km. Given that Gliese 
1214b has an M4.5 V primary, this implies that H/R^ ~ 0.0014. For an H to i?* ratio of this 
magnitude, ground-based measurements are presently capable of distinguishing between a 
clear if2-dominated solar-composition atmosphere, which will show measurable differences in 
planetary radius as a function of wavelength, A(i?p/i?*) ~ 0.005 (de Mooij et al. 2012), and 


a high-molecular weight atmosphere with a small scale height and a transmission spectrum 
that presents much smaller variations. 

Measurements of A(i?p/i?*) taken at different bandpasses during transits by Gliese 
1214b have indicated that the transmission spectrum is fairly flat between A = 0.6 /rm 


and A = 5 /xm ( 

Bean et al. 

2010 

Desert et al. 

2011; 

Groll et al. 

2011a 

Bean et al. 

2011 

Berta et al. 

2012 

de Mooij et al. 

2012 

). As of 2014, number of these measurements are 


in conflict, however, and while a clear solar-composition atmosphere is inconsistent with 
the data, Gliese 1214b could be presenting either a hazy low-metallicity H-He dominated 
atmosphere, or a low-scale height atmosphere composed of species (such as H2O) of high 


molecular weight. Kreidberg et al. (2014) report intensive HST observations of Gliese 1214b 
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(15 transits over 60 telescope orbits) that are sensitive enongh to reveal absorption featnres 
from constitnents snch as H2O, CO2, and NH4 in a high molecnlar weight atmosphere. A 
flat transmission spectrnm between 1.1 /xm and 1.7 /xm was observed, snggesting that Gliese 
1214b’s atmosphere must contain clonds or high-altitnde hazes. In addition, other planets 
with P < 100 d and M < 30 M0 are nearly certain to be discovered in transit across bright 


low-mass primaries. For example, Bonfils et al. (2012) have annonnced the detection of 
Gliese 3470b, which has Mp = 14 M0, Rp = 4.2 R0, and Tpi ~ 800 K, and which orbits an 
R* ~ 0.5 R© Ml.5 V primary. Spectrographic observations of of Gliese 3470b by Grossfield 


et al. (2013) nsing the Keck telescope are consistent with a flat spectrnm between 2.09 /xm 


and 2.36 /xm, which rnle ont a clond-free atmosphere in chemical eqnilibrinm, and which are 
again snggestive of high-altitnde haze or clonds. The prospects for improved characterization 
of planets in the nngiant category are especially bright, given that NASA’s forthcoming TESS 
Mission will locate effectively all of the transiting short period planets with radii Rp > 1.5 
R© orbiting stars with V < 12. 

Transit and occnltation measnrements act as effective bookends to full-phase photomet¬ 
ric measnrements, in which planetary emission is monitored throngh the conrse of an orbit. 


The most freqnently cited example of this type of observation was pnblished by Knntson 


et al. (2007), who presented a high-qnality 8/xm light cnrve for HD 189733b that spans the 


29-honr interval from jnst before the planet’s primary transit nntil jnst after its secondary 
eclipse. Adopting the assnmption that the planet is in spin-synchronons rotation with its 


orbit, Knntson et al. (2007) derived a global temperatnre distribntion for the planetary 


photosphere as a fnnction of longitnde, and fonnd that the hottest region of the planet is 
shifted east of the snbstellar point (see Gowan & Agol (2008) for the details of the map¬ 


ping techniqne). The longitndinal shift of the planet’s temperatnre maximnm is attribnted 
to snper-rotating v ~ lkms“^ eqnatorial winds, which are a robust feature that emerge 
from a wide variety of hydrodynamical models of the surface flows on tidally locked close- 
in extrasolar planets, (see, e.g. Showman et al. 2009; Showman & Polvani 2011[ ). In the 
past several years, a number of additional infrared phase curves in assorted bandpasses have 


been produced for a variety of strongly irradiated planets. Examples include Knutson et al. 


(2009b)’s follow-on partial-orbit MIPS 24/xm and Knutson et al. (2012)’s full-orbit 3.6/xm 


and 3.5/xm IRAG phase curves for HD 189733 (which corroborate the basic picture inferred 


by 

Knutson et al. 

2007 

), and a full-orbit optical phase curve from the Kepler satellite for 

HAT-P-7b (fBorucki et al. 

2009 

) that shows evidence of both dayside thermal emission and 


reflected light. Full or partial-orbit Near-IR phase curves have also been published for HD 


149026b (Knutson et al. 2009a), WASP-18b (Maxted et al. 2013), and WASP-12b (Gowan 


et al. 2012), all of which are hot, short-period planets on nearly circular orbits. The phase 


curves for these systems suggest that the selection of planets surveyed have a large range of 
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Time from Periapse (Earth days) 


Fig. 9.— Photometric IR phase curves for the eccentric planet HAT-P-2b, with orbital 


period, P=5.633d, eccentricity, e=0.511, and M = 9Mjup, adapted from Lewis et al. (2013). 
The top panel shows the geometry of the HAT-P-2b orbit, with periastron passage occurring 
14 hours after the primary transit. The lower panels display the light curves for the 3.6 
(top), 4.5/im (middle), and 8.0 /rm photometry (hlled circles). Optimized phase, transit. 


and secondary eclipse fits are over-plotted as a red line, see Lewis et al. (2013) for details. 
The dashed line represents the stellar flux level. 


efficiency in their ability to equilibrate day-night temperature distributions. 

A typical hot Jupiter with Teq ~ 1500K has its blackbody emission peak at A ~ 1.9/im, 
which is blueward of the Spitzer bandpasses. It is therefore useful to obtain phase curves 


near the region where a given planet disgorges the majority of its energy. Stevenson et al. 


(2014) used the HST WFC3 instrument to obtain full-phase 15-channel spectrophotometry 


in the l.l/im —1.7/im region for three consecutive full orbits of WASP 43b. This planet 
is a representative hot Jupiter, albeit with a cool (Teg = 4500, K7) primary, and a short 
P = 19.5h orbit. Stevenson et al.’s analysis indicates that the phase-resolved spectrum shows 
little contribution from molecular absorbers, with the only significant inffuence coming from 
H2O. In agreement with many of the other extant full-phase observations, the light curve 
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indicates that the planet’s hottest location is offset eastward from the substellar point, and 
that there is poor dayside - nightside heat circulation. 


Figure [^shows phase curves in several bandpasses for the eccentric planet HAT-P-2b, as 
presented by Lewis et ah (2013). Planets with eccentric orbits cannot rotate synchronously, 
and hence their atmospheres are virtually certain to display strongly time-dependent me¬ 
teorology. Phase-curve observations of such planets can potentially determine atmospheric 
radiative time constants, tidal luminosities, and pseudo-synchronous rotation rates, as well 
as the global manifestation of time-dependent storms (Laughlin et ah 2009). 


6. Planet Formation 

The vast majority of work on planetary formation has been done to explain the proper¬ 
ties of our Solar System, for which we have the most detailed data. We begin with a synopsis 
of these models, and then discuss some recent models of the formation and migration of ex¬ 
oplanets that travel on short-period orbits. 


6.1. An Overview Informed by the Solar System 


Models of Solar System formation are based upon data from our own Solar System as 
well as observations of the circumstellar environments of young stars. The nearly planar 
and almost circular orbits of the planets in our Solar System argue strongly for planetary 
formation within flattened circumstellar disks. When a molecular cloud core undergoes 
gravitational collapse, an occurrence triggered by loss of magnetic field and turbulent velocity 
support, its spin angular momentum leads to the formation of pressure-supported protostars 
surrounded by rotationally-supported disks, e.g., Terebey et ah (1984). Molecular cloud 


cores are sub-condensations of mass, M > Mq, characteristic size r ~ 0.25 pc, and rotation 
rate typically uj ~ 10“^^ Hz that are found within giant molecular cloud complexes in the 
galactic disk that typically contain 10"^ to 10® Mq. 


Circumstellar disks are commonly observed in association with newly formed stars (Ma- 


majek 2009), and are analogous to the primordial solar nebula that was initially conceived 


by Kant and Laplace to explain the observed properties of our Solar System; for a more 


detailed historical background, see Woolfson (1993). Planets form within such disks, and 


play a major role in the later stages of disk evolution. Terrestrial planets grow via pairwise 
accretion of solid bodies. In the generally favored core-accretion paradigm, giant planets 
begin their growth as do terrestrial planets, but they become massive enough to accumulate 
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substantial amounts of gas before the protoplanetary disk dissipates. Under some circum¬ 
stances, as described for example, by Gammie ( 2001 [), giant planets may form directly from 


the protostellar disk via gravitational instability (see, e.g. Durisen et ah 2007; Boley 2009), 
especially at large (d > 40AU) from the parent star. 

Within the Solar System, bodies up to the mass of Earth consist almost entirely of 
condensable (under reasonable protoplanetary disk conditions) material, and the fraction 
of highly volatile gas increases with mass above that of Earth. The most massive planets 
contain a considerable fraction of light gases. About 90% of Jupiter’s mass is H and He, 
and these two light elements make up ~ 75% of Saturn. In both cases, these planets contain 
~ 30 M 0 in elements other than hydrogen and helium. The large amounts of H and He 
contained in Jupiter and Saturn imply that these planets must have formed within ~ 10^ 
years of the collapse of the Solar System’s natal cloud, before the gas in the protoplanetary 
disk was swept away. The two largest planets in our Solar System are generally referred to 
as gas giants^ even though these elements aren’t gases at the high pressures that most of 
the material in Jupiter and Saturn is subjected to. Analogously, Uranus and Neptune are 
frequently referred to as ice giants, even though the astrophysical “ices” such as H 2 O, CH 4 , 
H 2 S and NH 3 that models suggest make up the majority of their mass ( Hubbard et al.||1995 ) 
are in fluid rather than solid form. Note that whereas H and He must make up the bulk 
of Jupiter and Saturn because no other elements can have such low densities at plausible 
temperatures, it is possible that Uranus and Neptune are primarily composed of a mixture 
of “rock” and H/He, although the high observed CH 4 /H 2 ratio in their atmospheres (30-40 x 
the solar value) makes this somewhat unlikely . This ambiguity in basic structure also holds 
for the majority of the ungiant exoplanets whose mass and radius have both been measured. 


Lithium and heavier elements constitute < 2% of the mass of a solar composition mix¬ 
ture. The atmospheric abundances of volatile gases heavier than helium (excluding neon, 
which is thought to be substantially depleted through gravitationally-induced settling) are 


3x solar in Jupiter (Young 2003), a bit more enriched in Saturn, and substantially more 


for Uranus and Neptune. The bulk enhancements in heavy elements relative to the solar 
value are roughly 5, 15, and 300 times for Jupiter, Saturn and Uranus/Neptune, respec¬ 
tively. Thus, all four giant planets accreted solid material substantially more effectively 
than gas from the surrounding nebula. Moreover, the total mass in heavy elements varies by 
only a factor of a few between the four planets, while the mass of H and He varies by about 
two orders of magnitude between Jupiter and Uranus/Neptune. 

There is a vast observational literature that details the presence, the properties, and the 


evolution of disks of Solar System dimensions around pre-main sequence stars (Williams & 


Cieza 2011). The existence of disks on scales of a few tens of astronomical units is inferred 
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from the power-law spectral energy distribution in the infrared over more than two orders 


of magnitude in wavelength (Chiang & Goldreich 1997 Haisch et ah 2001 Megeath et ah 


2004 Andrews et ah 2010). Observations of infrared excesses in the spectra of young stars 
suggest that the lifetimes of protoplanetary disks span the range of 10® - 10^ years (Strom 


et ah 1993 Alencar & Batalha 2002 Ercolano & Koepferl 2012). Much lower mass second 


generation debris disks composed of dust eroded from larger solid objects that accreted within 
protoplanetary disks are observed around main-sequence stars, with younger stars typically 


(but far from always) possessing more massive debris disks (Meyer et al. 2004 Krist et ah 


2012). Many debris disks appear to have inner holes, and lower mass stars typically retain 


their debris disks for more time (Wyatt 2008) 


Although to a hrst approximation the gas in the disk is centrifugally supported in balance 
with the star’s gravity, negative radial pressure gradients provide a small outwardly-directed 
force that acts to reduce the effective gravity, so the gas rotates about the star at slightly less 
than the keplerian velocity. Small solid bodies (dust grains) rotate with the gas. Large solid 
bodies orbit at the keplerian velocity, and medium-sized particles move at a rate intermediate 
between the gas velocity and the keplerian velocity; thus medium-sized and larger bodies are 
subjected to a headwind from the gas (Adachi et al. 1976| . This headwind removes angular 
momentum from the particles, causing them to spiral inwards towards the star. This inward 
drift can be very rapid, especially for particles whose coupling time to the gas is similar to 
their orbital period; smaller particles drift less rapidly because the headwind they face is 
not as strong, whereas large particles drift less because they have a greater mass to surface 
area ratio. Orbital decay times for meter-sized particles at 1 AU from the Sun have been 
estimated to be only ~100 years ( Weidenschilling||1977[ ). The large radial velocities of bodies 
in this size range relative to both larger and smaller particles implies frequent collisions, so 
most solid bodies may pass through the critical size range quickly without substantial radial 
drift. However, it is also possible that a large amount of solid planetary material is lost from 
the disk in this manner. Solid bodies larger than ~1 kilometer in size face a headwind only 
slightly faster than meter-size objects, and because of their much greater mass/surface area 
ratio they suffer far less orbital decay from interactions with the gas in their path. Thus, 
kilometer-sized or larger planetesimals appear to be reasonably safe from loss (unless they 
are ground down to smaller sizes via disruptive collisions) until some of these planetesimals 
grow into planetary-sized bodies. 


The initial composition of the protoplanetary disk is assumed to still be essentially 
identical to that of the protosun. Dust within a protoplanetary disk initially agglomerates 
via sticking/local electromagnetic forces that produce inelastic collisions. Particles gradually 
settle towards the disk mid-plane as they grow. 
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The conceptual bottleneck created by the short residence times of meter-sized particles 
in protoplanetary disks is perhaps the least well understood phase of solid planet formation. 
It is not yet clear how agglomeration from cm-sized pebbles to km-sized or larger bodies that 


are referred to as planetesimals occurs. Collective gravitational instabilities (Safronov 1969 


Goldreich & Ward 1973) might be important, although turbulence may prevent protoplane¬ 


tary dust layers from becoming thin enough to be gravitationally unstable (Weidenschilling 


& Cuzzi 1993). The interactions between small particles and gas in a rotating disk involve 


complicated fluid-particle interactions, and planetesimal formation is a very active research 
area. Recent work suggests that the answer may he in the joint action of aerodynamic 
and gravitational instabilities to agglomerate particles - see the review of Chiang fc Youdin 
(2010). Some models suggest that pebble-sized particles accumulated into planetesimals 


more than 100 km in size (Johansen et ah 2012) 


The primary perturbations on the keplerian orbits of kilometer-sized and larger bodies in 


protoplanetary disks are pairwise gravitational interactions and physical collisions (Safronov 


1969). These interactions lead to accretion (and in some cases erosion and fragmentation) of 
planetesimals. Gravitational encounters are able to stir planetesimal random velocities up 


to the escape speed from the largest common planetesimals in the swarm (Safronov 1969). 


The most massive planets have the largest gravitationally-enhanced collision cross-sections, 
and accrete almost everything with which they collide. If the random velocities of most 
planetesimals remain much smaller than the escape speed from the largest bodies, then 


these “planetary embryos” grow extremely rapidly (Greenzweig & Lissauer 1992). The size 


distribution of solid bodies becomes quite skewed, with a few large bodies growing much faster 


than the rest of the swarm in a process referred to as runaway accretion (Greenberg et ah 


1978 

Wetherill & Stewart 

1989 

). If embryos dominate the stirring, then M oc in this 

circumstance, if individua 

embryos control the velocity in their own zones, larger embryos 


take longer to double in mass than do smaller ones, although embryos of all masses continue 
their runaway growth relative to surrounding planetesimals; this phase of rapid accretion of 
planetary embryos is known as oligarchic growth (Kokubo fc Ida|1998). Eventually, planetary 


embryos accrete most of the (slowly moving) solids within their gravitational reach, and the 
runaway/oligarchic growth phase ends. 


The self-limiting nature of runaway and oligarchic growth implies that massive plan¬ 
etary embryos form at regular intervals in semi-major axis. The agglomeration of these 
embryos into a small number of widely spaced terrestrial planets necessarily requires a stage 
characterized by large orbital eccentricities, signihcant radial mixing, and giant impacts. At 
the end of the rapid-growth phase, most of the original mass is contained in the large bodies, 
so their random velocities are no longer strongly damped by energy equipartition with the 
smaller planetesimals. Mutual gravitational scattering can pump up the relative velocities 
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of the planetary embryos to values comparable to the surface escape velocity of the largest 
embryos, which is sufficient to ensure their mutual accumulation into planets. The large 
velocities imply small collision cross-sections and hence long accretion times. 


Once the planetary embryos have perturbed one another into crossing orbits, their 
subsequent orbital evolution is governed by close gravitational encounters and violent, highly 
inelastic collisions. Terrestrial planets continue to grow by pairwise accretion of solid bodies 
until the spacing of planetary orbits becomes sufficient for the conhguration to be stable to 


gravitational interactions among the planets for the lifetime of the system (Safronov 1969 


Lissauer 1995 Chambers|2001a Laskarj2000 Hansen &: Murray|2012). This process has been 


studied using iV-body integrations of planetary embryo orbits. Many of the simulations of 
this type have endeavored to reproduce our Solar System; they generally begin with about 2 
M 0 of material spread through the terrestrial planet zone, typically divided (not necessarily 
equally) among hundreds or thousands of bodies. These simulations generally display a high 
efficiency of converting seed embryos into terrestrial planets, and the end result is generally 
the formation of 2 - 5 terrestrial planets on a timescale of about 10^ years. Some of these 
systems look startlingly similar to our inner Solar System, with the caveat that the emergent 
planets often travel on more eccentric orbits that those of the terrestrial planets, although the 
simulations that begin with the largest numbers of bodies yield eccentricities more similar 
to those in our Solar System (O’Brien et al. 2006). It is possible that the Solar System is 
by chance near the quiescent end of the distribution of terrestrial planets, but more likely, a 
damping process of some type plays a role. Processes such as fragmentation and gravitational 
interactions with a remaining population of small debris, and interactions with a remnant 
gas disk (Ogihara et al. 2007) have been invoked to lower the characteristic eccentricities 
and inclinations of the ensemble of terrestrial planets. 


An important result of these iV-body simulations is that planetary embryo orbits execute 
a random walk in semi-major axis as a consequence of successive close encounters. The 
resulting widespread mixing of material throughout the terrestrial planet region diminishes 
any chemical gradients that may have existed when planetesimals formed, although some 
correlations between the final heliocentric distance of a planet and the region where most of 
its constituents originated are preserved in the simulations. Nonetheless, these dynamical 
studies imply that Mercury’s high iron abundance is unlikely to have arisen from chemical 
fractionation in the solar nebula. 


The mutual accumulation of numerous planetary embryos into a small number of plan¬ 
ets must have entailed many collisions between protoplanets of comparable size. Mercury’s 
original silicate mantle was probably partially stripped off by one or more of such giant im¬ 
pacts, leaving behind an iron-rich core. Accretion simulations also lend support to the giant 
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impact hypothesis for the origin of the Earth’s Moon; during the hnal stage of accumulation, 
an Earth-size planet is typically found to collide with several objects as large as the Moon 


and frequently one body as massive as Mars (Chambers 2001b). 


Models for the formation of gas giant planets were reviewed by Wuchterl et ah (2000) 


and Helled et ah (2013). Star-like direct quasi-spherical collapse is not considered viable, 


both because of the observed brown dwarf desert mentioned above and theoretical arguments 


against the formation of Jupiter-mass objects via fragmentation (Bodenheimer et ah 2000). 


The theory of giant planet formation that is favored by most researchers is the core-nucleated 
accretion model, in which the planet’s initial growth resembles that of a terrestrial planet, 
but eventually becomes sufficiently massive (several M 0 ) that it is able to accumulate sub¬ 
stantial amounts of gas from the surrounding protoplanetary disk. The only other hypothesis 
receiving significant attention is the gas instability model, in which the giant planet forms 
directly from the contraction of a clump that was produced via a gravitational instability in 


the protoplanetary disk (Durisen et ah 2007). Formation of giant planets via gas instability 


has never been demonstrated for realistic disk conditions. Moreover, this model has difficulty 
explaining the super-solar abundances of heavy elements in Jupiter and Saturn, and it does 
not explain the origin of planets like Uranus and Neptune. Nonetheless, it is possible that 
some giant planets within the currently observed census did form via disk instability. 


The core-nucleated accretion model (Mizuno 1980 Pollack et ah |1996), has been re¬ 


viewed by Lissauer & Stevenson (2007), with further details given in D’Angelo et ah (2011). 


The model relies on a combination of planetesimal accretion and gravitational accumulation 
of gas. In this theory, the core of the giant planet forms first by accretion of planetesimals, 
while only a small amount of gas is accreted. Core accretion rates depend upon the surface 


mass density of solids in the disk and physical assumptions regarding gas drag, etc. (Lissauer 


1987 Inaba et al. 2003). The escape velocity from a planetary embryo with M > 0.1 M 0 is 


larger than the sound speed in the gaseous protoplanetary disk. Such a growing planetary 
core hrst attains a quasi-static atmosphere that undergoes Kelvin-Helmholtz contraction 
as the energy released by the accretion of planetesimals and gas is radiated away at the 
photosphere. 

During the runaway planetesimal accretion epoch, the protoplanet’s mass increases 
rapidly. The internal temperature and thermal pressure increase as well, preventing sub¬ 
stantial amounts of nebular gas from falling onto the protoplanet. When the planetesimal 
accretion rate decreases, gas falls onto the protoplanet more rapidly. The protoplanet ac¬ 
cumulates gas at a gradually increasing rate until its gas mass is comparable to its heavy 
element mass. The key factor limiting gas accumulation during this phase of growth is the 
protoplanet’s ability to radiate away energy and contract (Pollack et al.||1996 Hubickyj et al. 
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2005). The rate of gas accretion then accelerates more rapidly, and a gas runaway occurs. 


The gas runaway continues as long as there is gas in the vicinity of the protoplanet’s orbit. 
The protoplanet may cut off its own supply of gas by gravitationally clearing a gap within the 
disk ( Lin fc Papaloizou||1979 ). Lissauer et ah (2009) used a three-dimensional adaptive mesh 
refinement code ( D’Angelo et al.||2002 2003) to follow the flow of gas onto an accreting giant 
planet, allowing the determination of final planetary mass as a function of the time-varying 
properties (density, temperature, viscosity, longevity, etc.) of the surrounding disk. Such a 
self-regulated growth limit provides a possible explanation to the observed mass distribution 
of extrasolar giant planets. 

Core-nucleated accretion also offers an explanation for the formation of low-density 
ungiant exoplanets. These bodies could be cores that grew massive enough to accrete some 
H 2 and He from their planetary disks, but were unable to accrete enough light gases to 


become giants (Rogers et ah 2011). 


6.2. How did short-period planets arrive at their observed locations? 


Quite aside from the zeroth-order questions related to their atmospheric properties and 
internal structures, the origins of the short-period planets are not fully explained. It is clear 
from Figure 2 that the hot Jupiters represent a population distinct from the short-period 
ungiants, and so it is reasonable to consider their formation scenarios separately. 

It is generally assumed that hot Jupiters were formed through the core accretion process, 
with assembly occurring at distances at least a few AU from their host stars, especially water 
ice, were available. As discussed above, in this standard core-accretion scenario for giant 
planet formation, a protoplanetary core composed largely of ices and rock reaches a mass of 
~ 10 M 0 , at which point it begins to rapidly accrete gas from the surrounding nebular disk. 
Thereafter, as a result of dynamical or hydrodynamical processes, the planet is delivered into 
a short period orbit. It is worth noting, however, that while disfavored, in situ formation 


for hot Jupiters has not been strictly ruled out. Detailed calculations by Bodenheimer et ah 


(2000), for example, suggest that hot Jupiters can form in place under certain conditions. 


The “classical” formation scenario for hot Jupiters was outlined by Lin et ah (1996), 


shortly after the discovery of 51 Peg b. In this picture, the planet grows to near its original 
mass at a distance several AU from the star. A planet with 51 Peg b’s mass is large enough 
to maintain a well-cleared annular gap in the disk that brackets the planet’s radial location. 


Once formed, the planet spirals inward via the process of Type II migration (Lin & Pa- 


paloizou 1986) on the parent protoplanetary disk’s viscous time scale, jv ~ 10® yr. 
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where v is the effective kinematic viscosity, and is the characteristic radial extent of the 
disk. The hnal location is often assnmed to be controlled either by disk truncation at an 
inner magnetospheric cavity, or via tidal interaction with the young, rapidly spinning parent 
star. This basic picture has been elaborated substantially during the past two decades. For 


overviews, see, e.g. Lubow & Ida (2011) and Kley & Nelson (2012). 


The disk migration scenario avoids the serious problems associated with in situ forma¬ 
tion, but is not itself without problems. One difficulty lies with so-called Type I migration 


Ward (1997), which is though to affect planets that have insufficient mass to clear and main¬ 
tain a gap in the protostellar disk. The timescale for Type I migration rj (Mp/M®)10® yr, 
is substantially shorter than the timescale that appears to be required to assemble a core 
large enough to initiate rapid gas accretion (Movshovitz et al. 2010). The mismatch hints 


at the possibility that a fundamental process is either missing or is not fully understood. A 
further problem arises from observed statistics of systems with massive, short-period planets. 


Hot Jupiters appear almost exclusively in single-planet configurations (Steffen et al. 2012) 
(with the V Andromedae planets providing a notable exception). 

Purely dynamical mechanisms have received increased attention in connection with the 


production of hot Jupiters. An early and influential study was carried out by (Rasio & Ford 


1996), who attributed hot Jupiters to the outcome of catastrophic planet-planet scattering. 


In their picture, extant short period planets are survivors of strong gravitational encounters, 
in which one planet was either ejected or placed in a long-period orbit, and the other is left 
in a high eccentric orbit with modest semi-major axis. Given a sufficiently close periastron 
distance, a high-eccentricity planet will subsequently undergo orbital circularization. 

The process of Lidov-Kozai migration with tidal friction has grown rapidly in popularity. 


In the Lidov-Kozai scenario, see e.g. (Eggleton & Kiseleva-Eggleton 2001 Wu & Murray 


2003 Fabrycky fc Tremaine||2007 ), a planet destined to become a hot Jupiter forms in a low- 
eccentricity orbit at a distance of several AU from a parent star of mass M*. In addition, the 
planetary orbit is substantially inclined at an initial angle io to the (for analytic purposes) 
circular orbit of a distant binary companion of mass M^. If the initial inclination angle, io, 
is larger (or better, substantially larger) than the Lidov-Kozai critical angle, io > Vnt = 
arccos[(3/5)^/^] = 39.2°, the planet is compelled to undergo Lidov-Kozai cycles in which 
the Jacobi energy, Ej = (1 — Cp)^/^cosip is conserved, and in which the planet experiences 
periodic cycling between successive states of high orbital eccentricity and high inclination to 
the binary plane. 


The associated timescale for these cycles is tlk = (2 T),^/(37rPp))((M* -|- Mp -|- Mb)/Mb). 
For large io, the planet’s periastron distance during the high-eccentricity phases of the cycle 
can be small enough so that dynamical tides acting on the planet efficiently convert orbital 
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energy into heat, thereby shrinking the semi-major axis, op, of the planetary orbit. The 
Lidov-Kozai cycles are eventually destroyed (or severely modihed) when the characteristic 
time scale for general relativistic precession, dictated by wgr = 3GM*np/(apc^(l — e|)) is 
substantially less than the Lidov-Kozai cycling timescale, i.e., 27r/a;GR < tlk- At this point 
the planet can transition to a long-running phase of secular evolution that is characterized 
by steady orbital decay and circularization, with the end product being a hot Jupiter. 


The Lidov-Kozai process is likely to deliver planets into orbits in which the planetary 
orbital angular momentum and the stellar spin angular momentum are initially largely uncor¬ 
related. A slew of measurements of the sky-projected angle. A, between the stellar rotation 
and planetary orbital angular momentum vectors using the Rossiter-McLaughlin spectro¬ 
scopic technique (Rossiter 1924 McLaughlin 1924) indicate that a substantial fraction of 


hot jupiter orbits are substantially misaligned with the equators of their parent stars (Winn 


& Fabrycky 2014), yielding evidence in favor of the Lidov-Kozai mechanism as a formation 


channel for a subset of the observed hot Jupiter population. 


7. Questions for the Next Ten Years and Beyond 

There are few precedents in the history of science in which a discipline moves so rapidly 
from shaky disrepute through a golden age of discovery and into a mature field of inquiry. In 
less than two decades, the study of extrasolar planets has accomplished all of this - answer¬ 
ing questions that were posed at the dawn of the scientific era, while affording tantalizing 
glimpses of revelations to come. This rapid progress partially obscures the fact that the ex¬ 
trasolar planets are fundamentally alien. Virtually none of their properties, either statistical 
or physical, aside from the scarcity of planets with 30 Mg < Mp < 100 Mg, were predicted 
or anticipated (Lissauer 1995, 2006), and theorists struggle to understand the zeroth-order 
features of the planetary distribution. At the current moment, it is thus useful to conclude 
by listing some provocative unsolved questions. 


What are the bulk compositions of the ungiants? If one considers logarithms of mass, 
the most prominent mass gaps among the bodies populating the Solar System are 
those between Jupiter and the Sun, with log]^o(Af©/Afj) = 3.02, and between Earth 
and Uranus, with logiQ^M^/M^) = 1.16. While it has long been possible to explore 
the properties of objects in the Sun-Jupiter gap through observations of low mass stars 
and brown dwarfs (Chabrier & Baraffe 2000), the mere presence of extrasolar planets 
with Mg < Mp < Mj orbiting main sequence stars has been known only since the 
discovery of Gliese 876 d, with Mp = 7.5 Mg ( Rivera et al.||2005 ), and the hrst density 
estimate became available only with the unveiling of Gliese 1214b, which has Mp = 6.3 
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Ma 


, Rp = 2.7 R 0 , and Pp = 1.7gem ^ Charbonneau et aT| (2009b). Structural models, 


as calculated by, e.g., Rogers & Seager (2010) or by Fortney et al. (2011), can easily 
reproduce Gliese 1214b’s mass and radius with configurations drawn from three board 
classes of models, (i) rock and iron planets with massive H/He or H 2 atmospheres, (ii) 
rock-ice cores surrounded by H-He envelopes of solar nebular composition, and (iii) 
so-called “water worlds” composed primarily of H 2 O. 


As we have described, attempts to characterize the atmosphere of Gliese 1214b have 
been met with ambiguity. It is likely that only when high-quality transmission spectra 
and interior mass distribution measurements are obtained that a clear understanding 
of the true range of compositions and structures of the ungiants will be obtained. 


What (if any) is the observable role of magnetohydrodynamics? Magnetic fields are 
important in nearly every branch of astrophysics. Earth and the four larger planets 
in our Solar System all harbor intrinsic magnetic fields that are generated by dynamo 
action in the planetary interiors. It is therefore natural to expect that many, perhaps 
most, short-period extrasolar planets with Mp > several M 0 also have magnetic fields, 
and to first approximation, these fields can be estimated by assuming that the Elsasser 
number, A = aiB"^/2pVt, is roughly constant across the aggregate of transiting planets 
ranging from Earth-like to Jupiter-like worlds. If we further assume synchronous ro¬ 
tation (which should hold for hot jupiters on circular orbits), then l/G cx The 

upper atmospheres of planets with P < 5 d are hot enough (T > 1200 K) that alkali 
metals (e.g., Na, K) have ionization fractions that are sufficient to couple the planetary 
magnetic field to the atmospheric hydrodynamic flow (see, e.g., Shu 1992). Studies, 
such as those by Batygin & Stevenson (2010), suggest that Ohmic dissipation is respon¬ 
sible for heating planetary interiors, thereby explaining much of the observed radius 
anomaly trend. Do MHD-mediated influences have a zeroth-order effect on exoplanet 
atmospheres, and can some of the bewildering variation in atmospheric emission be 
understood as arising from magnetohydrodynamic phenomena? 


• Where in the protostellar disks did the observed planets form? Our Solar System has 
long been thought to be the product of planet formation that occurred more or less 
in situ, via core accretion for the giant planets and via planetesimal accretion for the 
terrestrial planets. The detection of 51 Peg b, along with other early extrasolar planet 
discoveries, led to the suggestion, which soon became a paradigm, that long-distance 
planetary migration can, and does occur. At present, there are substantial arguments 
both in favor and against in situ formation of close-in ungiants, and for hot jupiters, 
where initial formation at large astrocentric distance seems somewhat more assured, 
it is not clear whether the LKGTF mechanism or quiescent long-distance Type-II disk 
migration is the dominant process. 
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Why is our inner Solar System empty? The Kepler Mission and the high-precision 
Doppler Velocity Snrveys have shown that the default mode of planet formation gen¬ 
erates multiple super-Earth mass planets with orbital periods, P < 100 d. Our own 
Solar System, however, contains nothing interior to Mercury’s 88-day orbit. Further¬ 
more, given the presence of Jupiter, with a mass and period combination that appear 
to be at least somewhat unusual, is there a possible connection between these two 


abnormal Solar System attributes? Work by (Walsh et ah 2011) suggests that both 
Jupiter and Saturn experienced significant radial migration, and the so-called “Nice 
Model” (Gomes et ah 2005; Morbidelli et ah 2005 Tsiganis et ah 2005) posits that 


substantial migration of the outer planets occurred during the solar system’s formation 
phases. An interesting, and as yet unanswered, question is whether such a process - 
generated, for example, via sweeping mean-motion resonances - could have thwarted 
planet formation in our inner Solar System, or destroyed planets that were originally 
there. 
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